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Abstract 

Searches for the permanent electric dipole moments (EDMs) of molecules, atoms, nucleons and 
nuclei provide powerful probes of CP violation both within and beyond the Standard Model (BSM). 
The interpretation of experimental EDM limits requires careful delineation of physics at a wide 
range of distance scales, from the long-range atomic and molecular scales to the short-distance 
dynamics of physics at or beyond the Fermi scale. In this review, we provide a framework for 
disentangling contributions from physics at these disparate scales, building out from the set of di- 
mension four and six effective operators that embody CP violation at the Fermi scale. We survey 
existing computations of hadronic and nuclear matrix elements associated with Fermi-scale CP 
violation in systems of experimental interest, and quantify the present level of theoretical uncer- 
tainty in these calculations. Using representative BSM scenarios of current interest, we illustrate 
how the interplay of physics at various scales generates EDMs at a potentially observable level. 
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1 Introduction 



Nuclear physics tests of fundamental symmetries have played a vital role in the development of the 
Standard Model (SM) and provide powerful probes of what may lie beyond it. As described elsewhere 
in this issue, these tests have uncovered the left-handed nature of the charged current weak interaction, 
helped single out the SM structure of the weak neutral current from various alternatives, revealed 
the phenomena of quark-mixing via the slight deviation from exact lepton-quark universality in weak 
decays, and provided stringent upper bounds on the neutrino mass scale. In this article, we focus on two 
symmetries for which nuclear physics studies have a long and illustrious history: time reversal invariance 
(T) and invariance under the combination of change conjugation (C) and parity (P). It is well-known, of 
course, that CP is not conserved in flavor- changing weak interactions, a phenomena now associated with 
the complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix. A consistent phenomenology 
of CKM CP violation (CPV) has emerged from extensive studies of K- and 5-meson properties and 
interactions. As a local quantum field theory satisfying the postulates of the CPT theorem, the SM 
thus also admits time-reversal violation, as the combined operation of CPT leaves the SM interactions 
unchanged. 

It is likely, however, that the SM picture of CP and T violation is incomplete. Numerous scenarios for 
physics beyond the SM (BSM) readily admit new sources of CPV. Given that the SM is likely embedded 
in a more complete theory of fundamental interactions, it is reasonable to expect novel signatures of CPV 
to appear along with other manifestations of new physics. Cosmology provides an additional compelling 
motivation for BSM CPV. Assuming the Universe was matter-antimatter symmetric at its birth or at 
the end of the inflationary epoch, additional sources of CPV are needed to explain the presently observed 
cosmic matter- antimatter asymmetry (for reviews and extensive references, see Refs. [IlElE])- From a 
perhaps even more speculative standpoint, the generation of the matter- antimatter asymmetry could 
also entail the violation of CPT invariance, while various exotic BSM frameworks also incorporate such 
a violation. While exploring the possibilities for BSM CPV and even CPT violation, one should bear in 
mind that there remains within the SM itself one as yet unobserved source of CP and T violation: the 
dimension-four QCD term, whose dimensionless coefficient, the vacuum angle 9, is now constrained 
to be no larger than ~ 10^^'^ by the non-observation of permanent electric dipole moments (EDMs) 
of the ^^^Hg atom and neutron. This exceedingly small upper limit - and the associated "strong CP 
problem" - has motivated the idea of an additional symmetry, the "Peccei-Quinn" (PQ) symmetry, 
whose spontaneous breakdown would imply the existence of the axion that has also not yet been 
observed. (For a recent review, see Ref. j4j.) 

In the quest to discover both BSM CPV as well as CPV generated by the SM strong interaction, 
EDM searches have generally provided by far the most powerful probes. In contrast to the CPV 
observed in the K- and i?-meson sectors, the existence of an EDM of an elementary particle or quantum 
system requires no flavor-changing interactions. The situation is more complicated in the SM, however, 
since CKM CPV requires the participation of three generations of quarks and, thus, flavor-changing 
interactions at the loop level. As a result, the EDMs of light quark and lepton systems generated by 
CKM CPV are highly suppressed. The individual quark EDMs vanish at two-loop order [5l E]. The 
lowest-order contribution to the neutron EDM, then, arises not from the individual quark EDMs but 
from a two-loop hadronic interaction involving two AS = 1 weak interactions between quarks. The CP- 
conserving interaction appears at tree-level, while the CPV interaction is generated from the one-loop 
d ^ s "penguin" operator that contains a sum over all three flavors of positive-charge quarks. The 
electron EDM first appears at four-loop order [7], suppressing it by several orders of magnitude with 
respect to the neutron EDM. 

Given the present and prospective EDM search sensitivities (see Table 22), one may consider CKM 
CPV to be something of a negligible "background", making these searches primarily probes of either 
SM strong or BSM CPV (together with P violation). Moreover, if the latter is fiavor-diagonal as one 
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encounters in many (but not all) BSM scenarios, then the sensitivity of EDMs generally exceeds that 
of other possible tests, such as CP-odd observables in a high-energy collider experiment. For these 
reasons, the emphasis in this article on nuclear physics tests of CP and T will fall on EDMs. 

Unraveling the implications of EDM searches is a multi-faceted problem, entailing physics at a 
variety of length scales. The experiments themselves are extraordinarily challenging, requiring exquisite 
control over a number of possible effects that could mimic an EDM. Still, significant improvements are 
expected for the traditional searches on neutral systems (Table [22| and proposals have been been made 
to measure the EDMs of charged particles in storage rings at similar levels. (For discussions of the 
experimental status and prospects, see, e.g. Ref. ^ as well as a forthcoming companion article to this 
review pj.) In what follows, we concentrate on the theoretical problem, seeking to provide a framework 
for interpreting experimental results that delineates the physics at different length scales that one must 
consider: the atomic, nuclear, and hadronic scales, wherein one contends with the complications of non- 
perturbative strong interactions and many-body physics; the Fermi scale, associated with the various 
effective operators outlined in the introductory article [lO]; the scale of BSM physics, 1/A, at which 
one encounters explicit new degrees of freedom whose interactions give rise to the effective-operator 
Wilson coefficients; the short-distance scale of high-energy collider experiments that may produce these 
new degrees of freedom directly; and the scales associated with early Universe cosmology that may be 
responsible for the generation of the matter-antimatter asymmetry. 

The interplay of these different scales is illustrated in Fig. [T] For purposes of this article, we assume 
the underlying dynamics of BSM CPV are associated with an energy scale A that lies well above the 
electroweak or Fermi scal^ If A < 10 TeV, high energy collider searches may discover the elementary 
particles responsible for BSM CPV, determine their masses, and provide information about the nature 
of their CP-conserving interactions. The new CPV interactions may also provide one of the ingredients 
needed for successful electroweak baryogenesis, though additional scalar degrees of freedom would also 
be expected in order to obtain a first order electroweak phase transition. In principle, collider searches 
could also observe the latter and measure their relevant properties [Ij. At the low-end of the energy 
scale, EDM searches look for the CPV "footprints" of these new interactions. In this energy regime, the 
extent to which the underlying CPV interactions become manifest depends on their interplay with the 
many-body and strong interaction dynamics of the hadronic, nuclear, atomic, and molecular systems of 
interest. The quantities that one extracts most directly from EDM searches, then, are not the underlying 
CPV interactions, but the hadronic, nuclear, and atomic matrix elements that they induce, such as the 
neutron EDM, time-reversal- violating and parity-violating (TVPV) nNN interaction, nuclear Schiff 
moment, and TVPV effective electron-nucleus interactiorj^ 

The bridge between these matrix elements and the underlying CPV dynamics is provided by a 
set of effective operators, whose coefficients are governed by an appropriate power of 1/A and a di- 
mensionless Wilson coefficient that depends on the details of the underlying dynamics. At energies 
between A and the weak scale, these operators contain all of the Standard Model fields and respect 
the SU(3)cxSU(2)2,xU(l)y gauge symmetry of the theory. Below the weak scale, the heavy SM de- 
grees of freedom are "integrated out", leaving a set of effective operators that respect the residual 
SU(3)cxU(1)em symmetry and that may have reduced mass dimension. Nevertheless, the Wilson co- 
efficients of the hadronic scale operators derive from those that enter at the weak scale, wherein the full 
gauge symmetry of the SM enforces certain relations between them. 

As we discuss below, it is in principle possible to use a combination of experimental results and 
theoretical hadronic, nuclear, and atomic computations to determine, or at least constrain, the Wilson 
coefficients without making any assumptions about the details of the underlying BSM CPV other than 

^It is possible that new CPV interactions are generated by new light degrees of freedom, a possibility that we do not 
treat extensively in this article. 

^Henceforth, we will use "CPV" when referring to the underlying elementary particle interactions and "TVPV" when 
referring to the resulting effects induced at the hadronic, nuclear, atomic, and molecular levels. 
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that it is associated with a high energy scale. In this sense, the theoretical effort associated with the 
bottom half of Fig. [T]is model- independent. More broadly, however, the goal of the "EDM program" is 
to derive as much information as possible from EDM searches, in conjunction with other precision tests 
and high-energy experiments, for both BSM and strong CPV as well as for the origin of matter. Doing 
so requires running operators from one scale to the next; matching the interactions at the boundaries 
between neighboring scales; and identifying and quantifying where possible the attendant theoretical 
uncertainties. One may then ask whether the emerging picture is consistent with any existing model 
for the underlying BSM CPV, precludes others, or perhaps points to one not yet invented. 

In what follows, we lay out the overall framework for this program. In doing so, we attempt to address 
a question that has been somewhat underemphasized in previous work, namely, the level of theoretical 
uncertainty associated with various steps in the interpretive chain of "running and matching". While 
the absence of theoretically robust hadronic, nuclear and atomic/molecular computations would not 
detract from the significance of the observation of an EDM, the level of theoretical uncertainty does 
affect one's ability to utilize present and prospective experimental results to pinpoint the underlying 
CPV mechanism or rule out various possibilities. In response to this issue, we provide a set of benchmark 
theoretical error bars associated with various quantities of interest, recognizing that this effort remains 
a work in progress and alerting the reader to the website where updated information will appear. 

The focus of this framework is on the ^-term as well as Wilson coefficients Ck for the dimension-six 
CPV operators indicated in the center of Fig. [Tj These operators, which break P as well as T, include the 
elementary fermion EDMs, the quark "chomo-electric dipole moments" (CEDMs), Weinberg three-gluon 
operator, and various four-fermion CPV operators (both semileptonic and non-leptonic). We provide 
a general set of expressions relating these operator coefficients to the hadronic, nuclear, atomic and 
molecular quantities of interest, given in Eqs. ( [3l2pl6| [4l68| [Kmpl73| and Table j?! We take into 



account the chirality-fiipping nature of the elementary fermion EDMs and quark CEDMs, writing the 
corresponding Wilson coefficients as the product of the fermion Yukawa couplings and a BSM scenario- 
dependent factor {5f or 5q). Doing so allows us to place the EDM and CEDM operator coefficients on the 
same footing as those for operators that do not flip chirality, such as the CPV electron-quark operators. 
A summary of the dimension-four and -six operators for light flavors (electron, up- and down-quarks, 
gluons) is given in Table ^ Note that one encounters thirteen quantities at this order, though some 
combinations have a more significant impact that others on the systems of experimental interest. One 
could, of course, expand the list to include the muon EDM as well as heavier flavors of quarks. In some 
cases the manifestation of the latter may in light quark systems may be non- negligible. Nonetheless, 
for purposes of this review we will concentrate largely on the already sizeable set of operators involving 
only the light flavor^ 

The dependence of various hadronic, atomic, and molecular quantities on 6, ImCfe, (5/, and 5q is 
then governed by the physics at the relevant scales. We compile the existing set of corresponding 
matrix-element calculations and give a set of benchmark values and theoretical ranges that can be used 
when extracting limits on 9, ImC^, 5f, and 5q from experimental results. For the discussion of hadronic 
matrix elements, we rely heavily on considerations of chiral symmetry as an overall guide, though we also 
quote results from lattice QCD, QCD sum rules, and quark models as well. An important conclusion 
from this survey is the need for a concerted future effort on the hadronic and nuclear matrix elements. 
While the literature on computations of the ImCfc, 5/, and 5q in various BSM scenarios is deep, the 
corresponding set of results for hadronic and nuclear matrix elements is relatively thin. Given the level 
of effort and resources devoted to the experimental measurements of EDMs, a commensurate attack on 
the theoretical side is clearly in order. 

Our discussion of this theoretical framework is organized in the remainder of the article as follows. In 

•^One should also bear in mind that not all possible sources of CPV naturally fit within the effective operator framework. 
If a new CPV interaction is mediated by a very light weakly coupled boson, the latter must be retained as an explicit 
degree of freedom 
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Figure 1: Electric dipole moments and the interplay of various scales. For purposes of 
illustration, only the impact of dimension six CPV operators is shown. Below the weak 
scale, some operators, such as the fermion EDMs and quark chromo EDMs are effectively 
dimension five, carrying an explicit factor of the Higgs vacuum expectation value {H^)- A 
summary of the operators of interest to this article appears in Table [TJ See text for a full 
discussion. 



Section [2| we briefly review the conventions and definitions, drawing on the notation of the introductory 
article. Section |3] contains a discussion of physics at the hadronic scale, including the running of the 
weak-scale operators to the hadronic scale, the various hadronic interactions cast in the context of 
chiral symmetry, and a summary of sensitivities of these hadronic quantities to the weak-scale operator 
coefficients. In Section [4], we review the status and open questions related to computations at the nuclear 
and atomic scales, including P- and T-odd nuclear moments such as the Schiff moment. We follow this 
discussion with an illustrative overview of the high-scale physics that may give rise to the weak-scale 
operators in Section [6j A discussion and outlook appears in Section [7j Throughout the article, we refer 
to other recent reviews [HI [12], [131 E] when appropriate, endeavoring to avoid excessively duplicating 
material that is amply covered elsewhere but updating when necessary. We also do not discuss other 
tests of CP and T violation, given the limitations of space for this review (for a discussion of T violation 
in neutron and nuclear /3-decay, see the companion article in this issue on charged current processes). 
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Table 1: Dimension four and dimension six CPV operator coefficients for light flavors. First column 
gives dimensionless Wilson coefficient (see Sec. |2]), followed by operator name and mass dimension 
(second column) and number of operators (third column). Final column indicates type of system in 
which a given operator will have its most significant impact. 



2 Conventions and Definitions 

The starting point for our analysis is the weak scale operators defined in the introductory article p^ . 
We concentrate on three sources of CPV, 

-^CPV = -^^CKM + Cfj + -C^Im • (2.1) 
Here the CPV SM CKM [iSj and QCD [ISl [IZl HE] interactions are 

£cKM = -^\/cMr^^I + h.c. , (2.2) 

C-e = -j|^^~Tr(G-^G,.) , (2.3) 

where g2 and are the weak and strong coupling constants, respectively, Uf^ (D^) is a generation-p 
left-handed up-type (down-type) quark field, Vq'^j^ denotes a CKM matrix element, Wj^ are the charged 
weak gauge fields, and G^j, = e^j,Q^G'^^/2 (eoi23 = 1 1^ is the dual to the gluon field strength G'^'^. In 
addition, 

^bIm = ^ (2-4) 

i 

gives the set of dimension-six CPV operators at the weak scale v = 246 GeV generated by BSM physics 
at a scale A> v. These operators \TU\ are listed in Tables |2] and |3} Note that the operators containing 
fermions are not CPV in and of themselves. Rather CPV effects arise when the corresponding coefficients 
Qfj-"''* are complex, as discussed below. 

''Note that our sign convention for e^i/Q.^, which foUows that of Ref. [T^, is opposite to what is used in Ref. [T3] and 
elsewhere. Consequently, Cg carries an overall —1 compared to what frequently appears in the literature. 
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Table 2: Dimension-six CPV operators involving gauge and/or Higgs degrees of freedom. Notation 
largely follows that of introductory article |10] with the following modifications: qp ^ Q denotes a 
left-handed quark doublet; F denotes a left-handed fermion doublet; denotes a right-handed SU(2) 
singlet fermion; and $ = (if) for being an up-type (down-type) fermion. For simplicity, generation 
indices have been omitted. 



{LR){RL) and {LR){LR) 


Qledq 


{UeR){dRQ^) 
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^lequ 
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Table 3: Dimension-six CPV operators involving four fermions. Notation as in Table |2j with the 
additional modification /p — )■ L with respect to the introductory article |10j . 



In this review we are mostly interested in the atomic/moleclar, hadronic, and nuclear aspects of 
CPV. We will therefore concentrate on the two lightest quarks, up and down, but will occasionally also 
point out effects of other quarks, especially strange. 



2.1 CPV at Dimension Four 



CPV from the ^-term in Eq. (2.3) is intimately connected with the quark masses. The "bar" notation 
indicates that this dimensionless quantity is a linear combination of a bare 6'-parameter and argument 
of the quark Yukawa couplings: 

^ = ^^ + argdet(A,Arf) , (2.5) 

where the second term arises after redefining the phases of all the quark fields. Alternatively, the 9- 
term can be eliminated, thanks to the axial anomaly, through a chiral rotation [20] ■ Enforcing vacuum 
stability to first order in the quark masses, all CPV is then in the quark bilinear 

Cq O £§pv = -m*e qi-i^q , (2.6) 

where ^ ^ 1 was used. Here 

m* = = — (1 - _ 2.7 

in terms of the average light quark mass m = {rriu + md)/2 and relative splitting e = (m^ — mu)/2m. 
Below we use = 2.3t.ol MeV and = 4.8tQ;3 MeV [21]. In the same notation the CP-even quark 
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mass operator is given by 



(2.8) 



Equation (2.6) will be the starting point in Section [3j for cliiral considerations that impact TVPV 
observables in nuclear physics. 



2.2 CPV at Dimension Six 

We now focus on a subset of the operators in Tables |2] and [3] that have been the the objects of most 
scrutiny, because they are expected to give the largest contributions at low energies [22]: Qgc (<? = u,d), 
Qfw, QfB, Qg, and various four-fermion operators. 

After electroweak symmetry breaking (EWSB) wherein cjp" — )■ (0, f /a/2), Qqc gives rise to the quark 
chromo-electric dipole moment (CEDM) interaction: 

^cEDM = ^ ^ ga-r^5g G;!, , (2.9) 

where = 1, • • • ,8) are the generators of the color group. Analogously, Q jw and QfB generate 

the elementary fermion EDM interactions, 

^EDM^_^^ (2.10) 

/ 



where F^y is the electromagnetic field strength. In the non-relativistic limit, Eq. (2.10) contains the 
CPV interactions with the electric field E, 

(2.11) 

/ 

where X/ is the Pauli spinor for fermion / and a is the vector of Pauli matrices. Thus, df gives the 
EDM typically quoted units of e cm or e fm. Letting 

a^^^guCfv,. (2.12) 

where Vk = B, W, and G for A; = 1,2,3 respectively, the relationships between the dq and df and the 
Cfv, are 

dq = ImC,G, (2.13) 

df = - 

where 



df = -V (xj I"^^/^' ^2.14) 



Im Cf^ = Im CfB + 2l[ Im Cfw , (2.15) 

and is the third component of weak isospin for fermion /. Here, we have expressed df and dq in terms 
of the Fermi scale 1/f , a dimensionless ratio involving the BSM scale A and v, and the dimensionless 
Wilson coefficients. Expressing these quantities in units of fm one has 

dq = -(1.13 X 10"^ fm) (J^y Im Cgc , (2.16) 
df = -(1.13 X 10"^ efm) (^^)^ Im C^^ . (2.17) 



It is useful to observe that the EDM and CEDM operator coefficients are typically proportional 
to the corresponding fermion masses, as the operators that generate them above the weak scale [Q^q, 
Qfw^ Q/b) contain explicit factors of the Higgs field dictated by electroweak gauge invariance. More 
physically, the EDM and CEDM operators - like the fermion magnetic moment - induce a flip of 
chirality and, thus, are naturally proportional to the fermion mass rnj = Yfv/\/2. Broadly speaking, 
then, one expects the Wilson coefficients to contain a factor of the fermion Yukawa coupling Yf. In the 
Minimial Supersymmetric Standard Model (MSSM), for example, a one-loop contribution to Im C^g 
from squark-gluino loop (see Fig. [s]) has the magnitude [23l [21] 

Im C,G = ^Y, sm[Arg{f,Msb*)] F{m,) , (2.18) 

where /i is the supersymmetric Higgs-Higgsino mass parameter, M3 is the soft SUSY-breaking gluino 
mass, b gives the Higgs SUSY-breaking mass parameter, and F{m,j) is a loop function that depends on 
the various superpartner masses rhj. In this case, the scale A would be the largest value of rhj entering 
the loop It is convenient to define two quantities Sq and Sf that embody all of the model-specific 
dynamics responsible for the EDM and CEDM apart from Yukawa insertion: 

ImC,G = Yq5q^dq = -—^ [£j 6q, (2.19) 

2 



2777 /■ / 1) ' 

Im Cf, ^ Yj6f^dj = -e-f (- j 5j . (2.20) 

While one often finds bounds on the elementary fermion EDM and CEDMs quoted in terms oi df and 
dq, the quantities 5/ and 5q are more appropriate when comparing with the Wilson coefficients of other 
dimension-six CPV operators, such as the three-gluon or semileptonic four-fermion interactions, that 
do not generally carry the Yukawa suppression. In what follows, we will provide expressions in terms 
of the ImC/y, (d/, dg), and (5/, 5q). In doing so, we will neglect the light-quark mass splitting and 
replace 

Y,-^Yq^^ (2.21) 

with ffi being the average light quark mass. 

The extraction of the CPV three-gluon and low-energy, fiavor-diagonal CPV four-fermion operators 
from is generally more straightforward. Making the identifications 

gives the so-called Weinberg three-gluon operator [25] 

nG _ 93 Cq rABCPiAiy^Bp^Cfi (c) 9q\ 



the CPV semileptonic interaction 



'^cpv = ^—^t^ [^756 dd - ee d-^^^d] - i^^^jY^ [eTse uu + ee (2.24) 

.(3) 



2A2 L J 



(^^uap ecr^ eua ''u , 



2A2 

^Note that we have not included a similar contribution involving the relative phases of A/3 and the squark triscalar 



terms (see Section 6.1) 
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and the CPV hadronic interaction [21] 



dd + uu d'y^d — d-y^u ud — du ^756?] 



+i 



2A2 

2A2^"'"^ [m75T^m dT^d + mT^m rf75T^rf - d-^^T'^u uT^d - dT^u u-y^T'^d] . 



(2.25) 



Note that in contrast to the all other CPV d = Q operators of interest here, the coefficient of the 
three-gluon operator (2.23) does not require the imaginary part. 



In addition to these four-fermion operators, the operator 

Qipud = i {(p^D^^if) UR'-f^dR 



(2.26) 



with Lf = ia24>* can also give rise to a four-fermion operator through exchange of the W boson contained 
in the covariant derivative. After EWSB, one has 



Q 



ipud 



(2.27) 



Exchange of the between the right-handed current in Eq. (2.27) and the left-handed current of the 

m 

(2.28) 



SM leads to an effective left-right (LR) Lagrangian with the CPV part given by [271 

^LR,CPV = -J^^^^^ [dLl^ULUR^^^dR - ULl^di dRy,jUR\ . 



After a Fierz transformation, one then obtains 



eff 

LR, CPV 



ImC^ud 
' 3A2 



[uu d'y^d — uy^u dd + 3 [uT^u dy^T^d — uy^T^u dT^d\ } 



(2.29) 



Although the RHS of Eq. (2.29) has the form of a product of scalar and pseudoscalar bilinears, it has a 
different flavor structure from the similar spacetime structures appearing in Eq. (2.25). As we discuss 
in Section [6| below, the interaction (2.29) is naturally generated in left-right symmetric theories. We 



(2.30) 



also note that the operator in Eq. (2.28) will mix [31] with an operator of the form 

dL-i^T^UL URj^T^dR - ULl^T^dL dR^^T^UR , 



generating the corresponding scalar ® pseudoscalar structures in Eq. (2.29), when running from the 
weak scale to the hadronic scale. 



2.3 Naturalness, Peccei-Quinn, and an Induced Vacuum Angle 



is^Hg EDMs imply that the coefficient of the 
In general, one would expect both terms of the 
right side of Eq. (2.5) to be 0{1). Obtaining a value that is ten or more orders of magnitude smaller 



It is well known that null results for the neutron and 
dimension four operator in Eq. (2.3) is tiny: 9 < lO"^''. 



would require a highly unnatural degree of fine-tuning to obtain a cancellation between the two terms. 



Note that in the limit of one massless quark, m* — )■ 0, CPV from Eq. (2.6) disappears entirely. However, 



such a possibility seems to be excluded on phenomenological grounds 

Alternatively, one may construct a mechanism that would generate a tiny ^ at a more fundamental 
level, through imposition of a symmetry or "geography" (see, e.g., Ref. [33]). The most well-known 
example of a symmetry argument is the Peccei-Quinn (PQ) mechanism. In brief, one starts from 
the anomalous axial U(l) symmetry of the SM in the limit of massless quarks, adding one or more 
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additional scalar fields to the SM whose interactions with the quarks preserves the tree-level axial U(l) 
symmetry. The enlarged symmetry, denoted U(1)pq, is spontaneously broken at a high scale, leading 
to a pseudoscalar Goldstone boson a, the axion. The corresponding axion Lagrangian is 



/:axion = ^d^ad.a - K(a) - ^ ^ Tr ( G^^G^, ) , (2.31) 



where the axion potential 



2 



n«) = 2X(0)(^^ + -^J +■■• (2.32) 

is given in terms of the topological susceptibihty x(0) as well as the axion decay decay constant fa 
whose value indicates the scale of spontaneous PQ-symmetry breaking. In two-flavor QCD, one finds 

x{0) = -m*{qq) , (2.33) 

with (qq) ^ — (225MeV)^ [H]. Physical observables depend on the combination 6 + (a)/ fa rather than 
on 6, where (a) is the axion vacuum expectation value (vev). Minimization of V{a) then implies that 
this combination vanishes, leading to a vanishing contribution to EDMs. The fluctuations about (a) 
correspond to the physical axion particle, whose mass is set by the ratio of \/xW) fa- 

Within the SM as well as in BSM scenarios, CPV radiative corrections to the quark masses (or 
Yukawa interactions) can generate a non- vanishing contribution to argdet(AttAd), re- introducing a pos- 
sibly unacceptably large magnitude for 6. If the given CPV scenario does not suppress these contribu- 
tions, the constraints on the underlying source of CPV can be quite severe. Invoking the PQ mechanism 
can alleviate these constraints. 

As emphasized in Ref. [H], the presence of higher dimension CPV operators Ocpv = ImC Qcpv/A^ 
can lead to an induced 6'-term. The operator Q^q in Table pi gives a tree-level shift in 9, which can 
still be removed through the PQ mechanism. More importantly, there is a shift in the minimum of the 
axion potential, which now reads [31] 



where 



Via) = x(0)oepv \^-jJ + 2^^^) V^Ja^ + ■ ■ ■ ' ^^-^^^ 



X(0)ocpv = lim I d^xe'""-" (0| T{GG{x), Ocpv(O)} |0) . (2.35) 
As a result the minimum of the potential occurs for 

^>^ = ^%^-^ina, (2.36) 

fa X(0) 

a so-called "induced" 6'-term|^ Thus, use of the PQ mechanism to eliminate the contribution of 9 to 
an EDM will introduce an additional contribution linear in the coefficient of a higher-dimensional CPV 
operator, ImC/A^. In the case of the CEDM operator, for example, one has 

^ind = ^Y.-^ (2-37) 
2 ^ mg 

q=u,d,s 

where ttIq characterizes the strength of the quark-gluon condensate {qo'^''T^G^^q) . In discussing the 
contributions of the dimension-six CPV operators to various P- and T-odd hadronic quantities, we will 
include the contributions from 6'ind wherever they have been explicitly computed. 



^Note that our definition gives an opposite sign to 6'ind compared to Ref. 
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3 CP and T at the Hadronic Scale 



In order to relate the interactions defined in Section [2] to P- and T-odd (TVPV)observables at the 
hadronic, nuclear, and atomic levels, we first introduce the most relevant hadronic quantities in the 
context of heavy baryon chiral perturbation theory (HBxPT) |35]. From the standpoint of effective 
field theory (EFT), HBxPT provides the natural and model-independent framework - consistent with 
the approximate chiral symmetry of QCD - in which to parameterize one's ignorance about presently in- 
calculable non-perturbative strong-interaction matrix elements of the various CPV operators appearing 
in /^cpv- For both the hadronic scale analysis as well as the Fermi scale effective operator formulation 



embodied in Eqs. (2.1), and (2.4), the EFT philosophy entails expressing the physical impact of un- 
known physics (BSM CPV or non-perturbative QCD) in terms of an infinite tower of operators having 
successively higher-mass dimension that carry appropriate inverse powers of the relevant mass scale: A 
in the case of BSM CPV and the QCD mass scale (or chiral symmetry-breaking scale) ~ 1 GeV 
in the case of the low-energy hadronic interaction. Doing so affords a systematic expansion of CPV 
observables in scale ratios, such as A^/A or -P/A^ where P denotes a soft momentum or pion mass. 
After truncation at a given order in these ratios, one has a reasonable estimate of the error incurred 
through omission of higher-order terms. 

Below all meson fields besides the pions can be "integrated out", their effects being captured 
by short-range interactions. Pions are light because they are the pseudo-Goldstone bosons of chiral 
symmetry, which plays an important role in determining the relative importance of the effective in- 
teractions. The term "heavy" in HB^PT indicates that one is only interested in dynamics where the 
nucleon is non-relativistic, having momentum = iriNV^ + k^ with being its velocity and [fc'^l <C ttitv, 
the nucleon mass. The nucleon is, then, described by a two-component field Ny{k) associated with a 
given velocity rather than a four-component Dirac field ipN- The anti- nucleon degrees of freedom are 
effectively also integrated out in terms of operators containing only Njj{k), its derivatives, and the pion 
field. Dropping the subscript "f " for notational simplicity, we give some representative terms in the 
resulting T- violating and P-violating (TVP V) Lagrangian [351 EZl [22] : 

ClY^ = -2iV (Jo + rfira) 5,Ar ^.F'^'^ + iV [#r ■ TT + ^7l^)7rV # (3r37r° ^ 

+CiNN [NS^'N) + C2NTN ■ df, {NS^tN) +■■■ . (3.38) 

Here, r and tt denote the isovectors of Pauli matrices and pion fields, respectively, while S'^ and f ^ 
denote the spin and velocity of the nucleon that take on values in the nucleon rest frame: S'^ — > (0, (?/2) 
when — )■ (1, 0). 



The first term in Eq. (3.38) defines the isoscalar (do) and isovector (di) "short-range" contributions 



to the nucleon EDM interaction: 

•HfoM = - {do ± Ji) X^(^X ■ E , (3.39) 

where the upper (lower) sign correspond to the proton (neutron) EDM interaction. 

The second term is the T- and P-odd pion- nucleon non-derivative interaction [38], consisting of 
isoscalar [g^^), isovector {glr^) and isotensor (f?!^^) pieces. These interactions have formally the same 
form when written in terms of a Dirac spinor tp^. Note, however, that various authors follow differing 
notation for the pion interactions. From the standpoint of a non-linear realization of chiral symmetry, it 
is more natural to build the Lagrangian from functions of tt/F^, where = 185 MeV is the pion decay 
constant. The resulting P- and T-odd ttNN couplings would then have dimension of mass. Moreover, 
in the absence of any breaking of chiral symmetry, the best choices of pion field are such that pion 
interactions involve derivatives, in which case the leading P-, T-even pion-nucleon interaction is of the 
pseudovector form. Yet, frequently a pseudoscalar form in terms of ipN is used. Chiral symmetry is then 
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(d^TT) ■ NrSm 


-2gA/F^ 


—gnNN/fnN 


—QtzNN /t^N 


-g/ruN 


-g/rriN 


-Iqa/F^ 



Table 4: Conventions for ttNN couplings: first three rows give TVPV non-derivative interactions while 
the last row gives the leading-order strong interaction. Note that denotes a heavy nucleon field, so 
the pseudovector interaction in the last row corresponds to the pseudoscalar coupling tt ■ ipNil^TipN 
(plus an additional two-pion interaction) in terms of a relativistic field ip^. In the chiral limit qa — 1-26 
and g = Qttnn — 12.6, while accounting for the Goldberger-Treiman discrepancy g = g-j^NN — 13.5 
[lallS] and gA ^ 1-33. Here, F^ = 185 MeV. 



only ensured if additional tt'^i^n'^Pn interactions are included. Fortunately, in most of the instances we 
are concerned with here, these additional interactions are irrelevant, and pseudoscalar and pseudovector 
interactions give the same result, once the corresponding couplings are related. A summary of notation 
used by various authors is given in Table |4} 

The third and fourth terms in Eq. ( |3.38[ ) contain T- and P-odd two-nucleon contact interactions, 
which represent all dynamics of range ~ such as vector meson (77, p, u, etc.) exchange. As 

we discuss below, these are expected to be the most significant short-range TVPV interactions among 
nucleons. The "■ ■ ■ " subsume an infinite number of other TVPV interactions: terms related to the 
above by chiral symmetry (see below) as well as interactions involving larger number of derivatives 
and nucleon fields and/or more powers of small parameters. For purposes of the present analysis we 
will not draw on these additional interactions explicitly. The reader should be warned that in general 
the Lagrangian (3.38) contains pion tadpoles [36], |22], as no spacetime symmetry forbids a tt" term 



(accompanied by its chiral partners with an odd number of pions) representing the disappearance of 
the neutral pion into vacuum. Tadpoles can be eliminated by field redefinitions, but for the left-right 
four-quark operator (2.28) a multi-pion vertex survives at leading order [22]. Although usually this is 



of no consequence, it might give rise to a significant TVPV three-nucleon force. 



The various hadronic interactions in Eq. (3.38) can be generated through the ^-term or any of the 



dimension-six CPV operators introduced above that contains only quarks and/or gluons. The semilep- 
tonic four-fermion operators Qeedq and Q^edu ^^^^ S^^^ rise to effective electron-hadron interactions. 
Concentrating on the electron-nucleon sector, we follow roughly the convention of Ref. [I3] to write 



r-eS 
'-eN 



V2 

-e^^a,3 ea^^e V^jycr"^ 



(1), 



CP + C'p^r, 



N 



(3.40) 



in terms of a relativistic nucleon field Normalizing to the Gp = ^/iV^v'^) allows us to make 

a straightforward comparison with limits quoted in the atomic EDM literature. A conversion to the 



operators normalized to A appears in Section 3.2.3 below. This interaction simplifies for a heavy nucleon 
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field, 



/-eff 
'-eN 



Gp r 

"71 



1 6^756 



N -8 ea^^e v^N 



a 



(0) 



+ ■ ■ 



(3.41) 



Here again, we neglect higher-derivative terms, operators containing more than two nucleon fields, and 
terms containing explicit factors of the pion field as implied by chiral symmetry. Note that the electron 
scalar nucleon pseudoscalar operators vanish at lowest order in the heavy baryon expansioiij^ 



3.1 Hadronic Matrix Elements 

In order to determine the dependence of the hadronic couplings defined above on the underlying sources 
of CP V, one must compute matrix elements of the various CP V operators introduced in Section |2j The 
result will be a set of expressions of the form 



d 



N 
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k 

k 



(3.42) 
(3.43) 
(3.44) 



where Ck denotes the Wilson coefficients for operator Qk-, as appropriate, and the coefficients etc. 
embody the results of the hadronic matrix-element computation. Note that for the three-gluon operator 
(2.23), here and in the rest of this review ImCfc stands for Cq. The coefficients and j3'P have the 



units of electric charge times length, and we will express all results as e fm. The coefficients A(j) and 
7^^'^^'' are dimensionless, while Ki and df'^ have dimensions of fm^. We note that presently very little is 
known about C7i_2. 

For future purposes, it will be convenient to define the sensitivity of the other hadronic quantities 



to either 5q and 5q or dq and dq via 



A 



V y 
A/ 



PflmCqG + PV^^Cq, 



75) ImC,G + 7?^Ima 
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e pj^ dq + Pjy dq 



A 



A 



(3.45) 
(3.46) 



Similarly, for the semileptonic interactions, we use Gp = 1/V2v'^; define g^^prp as the isoscalar and 
isovector form factors in the limit of isospsin symmetry 



{N\ [uTu + dTd] \N) 



(A^l [uTu - dVd] \N) = g^^^iPNTTsijN 



(3.47) 
(3.48) 



^Note that we have introduced an overall minus sign on the right hand sides of Eqs. ( 3.40|3.41 1 to match the convention 
in Ref. |13j and elsewhere, where the corresponding coefficients are defined for the Hamiltonian rather than the Lagrangian. 
Note also that an explicit —1 appears in front of the tensor interactions in order to facilitate comparison with other work 
in which an opposite sign convention is used for e^^ap- 
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where T = 1, 75, a^iy] and write for Cfp^rp, 

Cf= -Sf (I)' ImC<,-) and 4"= (^)^„,CW (3.49) 

4°'= g'? lmC!,:> and 4"= -9? Im4,-' (3.50) 



-9r(T)'lm41 and 4"= -si-" fy)' Im4f> . (3.51) 



where we define the combinations 

Cit^ = Cee,, ± C^ll . (3.52) 

For the dimension-six operators generated by BSM CPV, performing the hadronic computation 
entails two successive steps of running and matching. 

(i) One must first run the operators perturbatively from the BSM scale A to the weak scale. After 
integrating out the heavy SM degrees of freedom with appropriate matching, one then continues 
the running from the weak scale to the hadronic scale. The quantities Cqc, dq, 6, etc. are then 
defined at the hadronic scale A^^, ~ 1 GeV where nucleon matrix elements are then taken. They 
can be related to the quantities at the BSM scale A through an appropriate "i^-factor" , as in 

Im [g^Cqc] (AJ = K^g Im [gsCqc] (A) , ^^^(AJ = iT^cb^KA) , (3.53) 
[gs~S,]{A^) = KqclgsSqU), 

where we follow the convention in the literature and bundle the strong coupling with the Wilson 
coefficients ImC^ etc. The i^-factors then relate the product of and the Wilson coefficients at 
the two scales A and A^^. 

(ii) Second, one must compute the relevant matrix element at the hadronic scale utilizing non- 
perturbative methods. For the QCD 6'-term, only the second step is required. 

Carrying out the perturbative running is generally straightforward. In general, one must account 
for mixing among various operators. The full anomalous dimension matrix that takes into account the 
EDM, CEDM, three-gluon and four quark operators has recently been obtained in Ref. [35] • Prior to this 
work, efforts concentrated largely on the evolution of the EDM, CEDM, and three-gluon operators P6| 
Wl\ HHl HH]. Within this limited subset of operators, only the three-gluon operator is multiplicatively 
renormalized. The resulting "ii"-factor" , obtained after taking into account two-loop running and 
threshold corrections, is given in the first line of Table [5j The three-gluon operator, however, will mix 
into the CEDM while the latter will mix into the EDM. Consequently it is not generally possible to 
quote a single A'-factor for the latter two operators. Since the work of Ref. |18], however, it has often 
been the practice to do so in the literature. The reason is that that in the MSSM, Cq arises at two-loop 
order, whereas the CEDM first occurs at one-loop. Thus, the mixing of Qq into Qgc is effectively higher 
loop order. In the more general case, one must consider the full effects of operator mixing. Nevertheless, 
for illustrative purposes we quote a i^-factor for the CEDM to illustrate the magnitude impact made by 
QCD evolution from the weak to hadronic scales. Under similar assumptions, the authors of Ref. |18] 
obtained the i^"-factor given in Table |5} 

Performing non-perturbative computations is a more challenging task. Before reviewing the status 
of such calculations, it is useful to delineate expectations for the hadronic matrix elements based on the 
chiral symmetry properties of the operators, following the framework developed in Refs. [361 122] . For 
a parallel treatment in the context of chiral SU(3), see Refs. [50] [51]. We follow this discussion with a 
review of explicit computations utilizing various approaches. 
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Operator 




Reference 


Remarks 


Qn 
^ G 


3.30 


Ii6l liTl H8l HQ] 


Mult Renorm 




3.30 


iS] 

! ! 


Mixing neglected 


Q,v,V = B,W 


1.53 


HE] 


Mixing neglected 


quqd 


matrix 







Table 5: Illustrative perturbative renormalization factors for dimension-six CPV operators, Q. As 
in Eq. (3.53) the i^-factors apply to the product of with the ImC^, etc. In general, only Qq is 
multiplicatively renormalized, with the renormalization if-factor given in the first row. For all other 
operators, one must take into account mixing. Under the assumptions made in Ref. |48j for the MSSM, 
approximate /^-factor for the EDM and CEDM operators have been obtained as quoted above. No 
analogous approximation has been made for the four-quark operators, so we do not list corresponding 
entries. For a recent determination of the anomalous dimension matrix, see Ref. ^5] . 



3.2 Chiral Symmetry and Naive Dimensional Analysis 



In the limit of vanishing quark masses, the QCD Lagrangian is invariant under separate rotations of 
the right- and left-handed fields. Specializing to the two lightest flavors, these rotations are given by 



q exp [it ■ {OrPr + OlPl)] q 



(3.54) 



where Pr(l) denote right- (left-) handed projection operators and 0b.{l) are three-component vectors of 



arbitrary real numbers. For future reference, it is useful to re-express Eq. (3.54) in terms of vector and 
axial rotations: 

q exp [it ■ {6v + OaI^)] q ■ (3.55) 



The chiral SU(2)/jxSU(2)2, transformation embodied in Eqs. (3.54, 3.55) are isomorphic to those of 
SO (4), and for present purposes it is convenient to consider objects that have definite SO (4) transfor- 
mation properties. For example, four-component SO (4) vectors 



V 



V 



change, under an infinitesimal transform, by 



sv 



V xOv + V^Oa 
-OaV 



(3.56) 



(3.57) 



where 6v,a are presumed to be tiny. 

Terms in the effective Lagrangian just above the hadronic scale, and in particular the CPV inter- 



actions in Eq. (2.1), break chiral symmetry in specific ways. In order to reproduce the corresponding 



S matrix, the effective Lagrangian written in terms of hadronic fields has interactions that break the 
symmetries in the same way. For example, instead of a component of a chiral four- vector V[q] built out 
of quark fields, there will be a hadronic chiral four- vector V^[7r, A^] built from nucleon and pion fields. 
For a particular choice of pion fields, the latter can be related to one having no pions, \^[0, A^], by 



V[7r,A^] 



V[Q,N]- 
V4[0,Ar]- 



277 



77 



DF^ V P 

277 / 77 



DF^ \F^ 



-Va%N] + V[^,N] 



(3.58) 
(3.59) 
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where D = 1 + ■n'^/F^. 

The proportionahty constant between the hadronic interaction strength and the interaction strength 
above the hadronic scale is the hadronic matrix element. Hadronic interactions obtained from different 
components of the same object share the same matrix element. When the matrix element is not 
known, it can be estimated using naive dimensional analysis (NDA) [521 [25] . Because the short- distance 
physics incorporated in operators in the Lagrangian cannot be separated from quantum-mechanical 
effects represented by loops in Feynman diagrams, one assumes that the natural size of the operator 
coefficients is given by loop cutoff changes of 0{1). If M denotes the scale of breakdown of the EFT, the 
dimensionless "reduced" coefficient [A-k)"^'^ M^~^g of an operator of canonical dimension D involving 
N fields is assumed to be C(l) times the appropriate powers of the reduced couplings of the underlying 
theory. When applied to chiral-symmetric operators, which are characterized only by the reduced 
QCD coupling gs/An, but to any power, consistency requires that we take ~ Att in matrix-element 
estimates. 

Before proceeding with detailed applications, we consider two simple illustrations. First, when 
matching the CEDM operators onto the nucleon EDMs, we note that both operators posses the same 
canonical dimension. In this case = 3 and D = 5 (as we are below the weak scale). Since dq and 
have dimension M~^, we need only focus on the powers g^ and Air. In this case, = 3 implies that 

pr^ = W47r~l . (3.60) 

On the other hand for the contribution of the three-gluon operator, which is D = 6 even below the 
weak scale, we require one additional power of the hadronic scale as in 

= g,Aj4n ~ A^. (3.61) 
A second simple example is provided by the quark mass effects on the proton and neutron masses. 



The quark mass operators in Eq. (2.8) are components of two S0(4) vectors whose fourth component 



transform scalar or pseudoscalar under parity. 



S[,] = { and P[g] = .^7^ ) , (3.62) 



respectively. Replacing the light quark doublet in Eq. (3.62) by heavy nucleon fields and noting that the 
pseudoscalar operators vanish to lowest order in the heavy nucleon limit, we obtain the corresponding 
nucleon SO (4) vectors: 

^[0, N] = ( ) and P[0, N] = (^l^\ . (3.63) 



Using Eqs. (3.58, 3.59) one finds that the fourth component of S[k, N] and third component of P[k, N] 
give contributions from the average quark mass and mass splitting to the average nucleon mass, {'ffiN)qi 
and nucleon mass difference, (Am7v)q = {rUn — rnp)q, respectively: 

^Iss = - {mN)q NN + '-^^Nt,N . (3.64) 

The reduced coefficients are {mj^)q/A^ and {Amiy)^/ A^, which should be linear in rh/A^ and erh/A^, 
respectively, so that from NDA one expects {ffii^)q ~ rfi and {Amj^)q ~ efh. 

These terms are linked by chiral symmetry to others that contain an even number of pion fields, which 
contribute to pionic processes such as pion-nucleon scattering and pion production in nucleon-nucleon 
collisions. The two terms in Eq. [3.64 can be seen as the isospin-symmetric and breaking components of 
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the sigma term. The corresponding coefficients have therefore been evaluated in lattice QCD and also 
extracted from data. Results are, by and large, in agreement. The extrapolation of lattice results on 
octet baryon masses in 2+1 flavor QCD [53l [M], for example, gives {jnN)q = 45 ± 6 MeV, which agrees 
with the venerable value from Ref. [55]. Other extractions from data give similar values (see, e.g., the 
compilation in Ref. [56]; see also [57] ). Similarly, the lattice value (AmAr)^ = 2.26 ± 0.57 ± 0.42 ± 0.10 
MeV [58] is consistent with other lattice evaluations (see Ref. [5^), with a determination of the 
electromagnetic splitting using dispersion relations |60], and with an extraction from pion production 

For future purposes, it is useful to relate m to the pion mass. Since the pion is the pseudo-Goldstone 
boson of spontaneously broken chiral symmetry, its mass vanishes in the limit m — )■ 0. Away from this 



limit, the pion mass term arises from a scalar SO (4) vector as in Eq. (3.63), but with A^A^ — )■ 1. It is 
thus proportional to m, and from NDA ~ mA-^. It is convenient, as we do below, to estimate the 
coefficients of chiral-breaking operators steming from the quark masses in terms of m^/A^. 

3.2.1 Applications: ^-term 

We now use these classifications to identify the expected scaling of various hadronic operators as they 



are generated by underlying CPV interactions. We begin with the QCD 0-term, Cq (2.3), in the form 



of the quark bilinear (2.6) 



This transformed 6'-term is the fourth component of the same chiral SO (4) pseudovector as the 



quark mass term, and it gives rise to the —2Nt ■ tzN/F.^ term in Eq. (3.38) [36j. Thus, Nt^N 
and —2Nt ■ ttN/Ft^ (plus its chiral partners with more pion fields) transform as the third and fourth 
components of the same SO (4) pseudovector at the hadronic level, and the coefficient of the latter must 
be given in terms of matrix elements of the former: 

ChadP[7r,A^] = (had|P|had) , (3.65) 

where the state |had) contains appropriate nucleon and pion modes, viz. 

(ATlPglAT) = ChadiVrgAT and {N\P^\N'k) = -2C^^Nt ■ -kN / . (3.66) 

On the other hand, since 

em{N\P,\N) = ^ '" Nt.N , (3.67) 

we have that 

Chad = (Am;v)g /2em . (3.68) 
Hence, the matrix element of the QCD ^-term operator is 

- ^ (1 - e^) ^ {N\P^\Nti) = Qj^^ . = ^WiVr ■ tzN . (3.69) 

2 2e F^ 

Thus, we obtain the prediction 

3V = ^^—(^ or A(o) = — — (3.70) 

If the matrix element (1 — e^)(AmAr)g/e is calculated, g'^'^ /9 comes for free. As a rough estimate, taking 
{AmN)g ~ 2em and fh ~ -F7r/20 one would expect gl^^ ~ 0.05^, or A(o) ~ 0.05. We may also express 
the relationships in Eq. ( 3.70[ ) as 

A,o, ~ ^ . (3.71) 
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This expectation is given in Table [15] along with predictions for the same quantity using other ap- 



proaches. The lattice value for (Amjv)g given above implies 

A(o) = 0.017 ±0.005 , (3.72) 

where the error is obtained by adding the lattice uncertainties in (AmAr)g in quadrature. 

The foregoing reasoning leads [36] to analogous expectations for the other A(j) as well as the hadronic 

coefficients oat, and 'y'^^'^y For example, the simplest way to produce the isovector TVPV ttNN 
interaction n^NN in Eq. (3.38) is from a tensor product of two pseudoscalar vectors, and as a conse- 
quence 

4 

TTl 

A(i) - ^ , (3.73) 



where we took e ~ 1. The isotensor ttNN interaction in Eq. (3.38) is even more suppressed. 

The analogous arguments for the short-range components of the nucleon EDM are more complicated 
because one needs to account for the chiral transformation properties of the interaction between quarks 
and the photon field A^, 

C\tc = -^^.?-7^(l + 3r3)g. (3.74) 

While the first term is a chiral scalar, the second is the 3-4 component of an antisymmetric tensor. 
Taking the tensor product with the pseudoscalar vector P, they give rise, respectively, to the isoscalar 
and isovector nucleon EDMs. Thus, one expects 

2 2 2 

4i~e^Tf oi a^~e-^~0.2-f efm, (3.75) 

X X X 

where the additional factors of are simply a consequence of dimensional analysis. 
3.2.2 Dimension-six operators 



We now consider the dimension-six CPV operators appearing in Eq. (2.4) and arising from BSM 



physics [22]. The quark CEDMs can be embedded in S0(4) vectors and pseudovectors: 

Thus, the isoscalar and isovector CEDM operators transform as the P4 and 5*3 components of an SO (4) 
pseudovector and vector, respectively. They contribute to g^^ and g^^ without any additional factors 
associated with chiral symmetry breaking. Thus, we expect these two couplings to be comparable. 



(2) 

Again (jj^ is a. higher-order effect. However, since the EDM and CEDM Wilson coefficients carry an 
explicit factor of the quark Yukawa coupli ngs Yg ~ rrig/v ~ m^/(f A^), it is useful to express the glr^ in 
terms of the quantity 6q appearing in Eq. (2.19) as well as the pion mass and QCD mass scales: 



-(0.1) ^l^X r ~q ^l^X (o 70\ 

~ 1^ ^(0,1) ~ ■ (3-78) 



Similarly, we obtain for the sensitivity of (^l^'^'* to the d, 
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From similar considerations, the CEDM contributions to the nucleon EDM requires no additional 
chiral suppression but only an electromagnetic interaction, leading to 



or, alternatively. 



eml ~ - ml 



<<(oa,~^*, or Co,,~^ (3,81) 



and 



~p% = l. (3.82) 
For the effect of the quark EDM operators, the logic is similar. The transformation properties of the 



isoscalar and isovector quark EDMs are obtained by replacing the T"^ G"^^ in Eq. (3.76 ) by F^'^ . Purely 



hadronic operators now require integrating out high-momentum photons exchanged among quarks, 
which generates an additional factor of at least a/vr. In particular, for the g^^ we expect 



~ ^ , 0Ju\ ~ ^ , and rif-^ ~ — - — ^ . (3.83) 



This suppression renders these operators irrelevant for most purposes. In contrast, no such a/vr factor 
is needed for the quark EDM contribution to the nucleon EDMs since the photon is external. In this 
case, we expect 

I^N ~ P'n ' Pn^Pn^ and C(o,i) ~ C(o,i) • (3.84) 

(1 8) 

The situation for the three-gluon and four-quark operators operators is more subtle. Both the Q^;'°^ 
and Qq are chiral-invariant pseudoscalars. Using S and P from Eq. (3.62) we obtain 



S ■ P = -qri-f^q ■ qrq + qq qi-f^q = QJJ^^^^ . (3.85) 

A similar result applies to the analogous definitions of S and P but with the SU(3)c generators 

(s) 

included, yielding the operator Q^J^^. The CPV three-gluon operator is trivially a chiral pseudoscalar 
as the gluon fields do not transform under SU(2)rxSU(2)l. Consequently, the contributions from these 
operators to the cjn^ require explicit factors of m^/A^ and em^/A^ that reflect the chiral symmetry 
breaking needed to generate components of an SO (4) vector. Letting Ck denote the Wilson coefficient 
for any one of these three operators we then expect 

or ,-^(^)^|, (3,86) 



_(2) 

with (yf^r yet again at a higher order. Note that for such chiral-invariant TVPV sources the contact 
interactions in Eq. (3.38), which are chiral invariant, can be generated without any suppression from 
ffi: 

''u~;P^ImC. or S^l ~ ^ ■ (3,87) 
In nuclei they can be competitive with one-pion exchange with one glr''^^ interaction, as the enhancement 



~ m~ from a pion propagator is compensated by the in Eq. (3.86). 

There is no suppression also when, in combination with the quark electromagnetic interaction, these 
sources produce isoscalar and isovector components of the nucleon EDM transforming as chiral scalar 
and antisymmetric tensor, respectively: 

rf(o,i)~^ImC, or (3-88) 
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Matching the four-quark operator in Eq. (2.29) onto hadronic operators, we first observe that it 
is the 3-4 component of a symmetric chiral tensor. Looking at the terms without color matrices T^, 
one may decompose them into two terms, each having the form of the product of the third and fourth 
components of either S or P: 

uu di'j^d — ui'y^u dd = - [S3 ® 5*4 + P4 (g) P3] . (3.89) 

Similar relations hold for terms with T^. 

The 5*3 ® 5*4 structure leads directly to a contribution to gl^\ However, the presence of signals 
vacuum instability, as it generates also terms with an odd number of pions, including a neutral pion 
tadpole. Such a tadpole can be eliminated by a chiral rotation, but it leaves behind a three-pion 
interaction that has not been well studied and is neglected here. After this rotation, one finds also 
contributions to both glr'^ and and 



As for CEDM, there is no suppression in the nucleon EDM 



- eA^ ImCy^rf >^ud) e^x /oqi^ 

3.2.3 Semileptonic interactions 

For the semileptonic matrix elements in Eqs. (3.49|3.50|3.51 ), we can follow similar considerations. 

'''''' / -\ 

Scalar and tensor interactions are most important at low energies, Eq. (3.41). The values of follow 
straightforwardly from the contributions of the light quarks to the nucleon masses (3.64): 

(0) ^ (m^ (1) ^ (A!^ . (3.92) 

2m ' Aerfi ^ ' 



From the empirical and lattice values for {m]^\ and (AmAr)g given in Section 3.2, one obtains 

^5°^ = 6.3 ± 0.8, ^^^^ = 0.45 ±0.15 , (3.93) 

where the errors do not include the range of values for the light quark masses. 
Similarly, the tensor matrix elements are related to those of the quark EDM, 

gf = Ap\,^. #=4p?,). (3.94) 

In principle, we could adopt the chiral SO (4) approach to estimate the g"p . Instead, we find it useful 
to follow a variant of the arguments discussed in Refs. [121 E2]- Taking the divergence of 

(A^|m7/,75M - di'^^'^^d\N) = gAi'Nlfilbni'N , (3.95) 

employing partial conservation of the isovector axial current, using nucleon equations of motion, and 
observing that the d ivergence of the isovector axial current is anomaly-free as well as the definitions of 
the gp^ in Eq. ( |3.47 ) we obtain 

9? = ^ , (3.96) 
m 

assuming exact isospin symmetry. Analogously, we can define for matrix elements of the remaining 
flavors 

{N\QTQ\N)=g^^NT^PN , (3.97) 
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with Q denoting here s, c, b, or t. We will henceforth treat the strange quark separately and refer to c, 
b, and t as the heavy flavors. 

The scalar form factors can be obtained from analyses of the a term. Defining 



{N\msss \N) = asNN 



and 



220 MeV 

as has been conventional in the recent literature (see, e.g. [I^), one has 



9s ~ 



220 MeV 



(3.98) 
(3.99) 

(3.100) 



with Ks = 21 ±6 MeV |5l]- For heavier flavors a short derivation appears in Appendix |A] and we simply 
quote the result here: 

o /66MeV\ , , , , 

g'^= 1-0.25/s:, . 3.101 

\ J 

Note that the strange and heavy quark scalar form factors are isoscalar. 

Again, the derivation for the pseudoscalar form factors is given in Appendix |A] and leads to 



{N\mQQi-f5Q\N) =ijN 



Q(0) Q(l) 

9p + 9p ^3 



N 



with 



9p 



9^ 



(0) 

9 a 



79a 



2mr 



+ 9a 



m,, 



2m, 



Q 



(3.102) 

(3.103) 
(3.104) 



With these considerations in mind, we now review explicit computations of the hadronic matrix 
elements, referring to the expectations based on considerations of chiral symmetry wherever possible. 
To date, first-principles computations have been undertaken for the quantity a„; the scalar coupling 
g^g'^^ that can be obtained from (Tttat and the contribution from the quark mass difference to the nucleon 
mass splitting as discussed above; direct computations of the isovector form factors g^"* for V = S,T 
using lattice method^for a compilation and recent results, see Ref. [HS! ); and for the nucleon electric 
dipole moment form factors using HBxPT. Indeed, HBxPT can provide considerable insight into the 
dependence of matrix elements on light quark masses while implementing a consistent expansion of QCD 
in scale ratios as discussed above. On the other hand, knowledge of the low-energy constants requires 
additional input, either from a lattice computation, direct measurement, or a model estimate. The 
present state of the art still relies heavily on the latter approach. As a result, there exists a considerable 
degree of model-dependent uncertainty in the values of the j3^^\ etc.. One objective of this review is to 
provide a set of benchmark values and theoretical uncertainties for these parameters. A compilation of 



existing results drawn from various methods is given in Tables 15 through 20 appearing in Appendix 
[Bj We now review explicit computations of the hadronic matrix elements, leading to the results quoted 
in these tables. 



^Note that the computations in Ref. 63J apply to the charged current form factors that are related to those of interest 
here by an isospin rotation. Note also that the relative normalizations of the quantities here and in that work are given 
by 3S,T = 2g^^^ . 
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3.3 Chiral Perturbation Theory 

At low energies QCD reduces to HBxPT, where the symmetries of QCD, including the considerations 



of Section 3.2, are naturally expressed. The isoscalar and isovector electric dipole form factors (EDFF) 
of the nucleon are defined through the TVPV electric current 



S^'vq-S- qv^ + — {S^'k -q-S- qk") + ■■■ 



where the outgoing photon momentum q = p — p', and k = {p + p')/^. The EDFF 

F,{Q^) = d,-S[Q' + --- 



(3.105) 



(3.106) 



gives the EDM at zero momentum transfer, = 0, and the term linear in (the form-factor radius) 
provides an electromagnetic contribution that cannot be separated from a short-range electron-nucleon 
interaction. 

HB^PT provides the momentum and pion mass dependence of the EDFF in terms of the pa- 



rameters appearing in the TVPV Lagrangian (3.38) and in the Lagrangian encoding P-even, T-even 



interactions. The full nucleon EDFF has been calculated to leading order (LO) [Ml ESI |66l [671 [22] and 
next-to-leading order (NLO) [681 [Ml [SZl UBl [22] in the P/Mqqd expansion for all sources described 
above: the ^-term and the dimension-six CPV operators. (For related results in three-flavor xPT, see 
Refs. [m [721 [721 IZH I5H1 EI]-) The resulting contributions fall into two classes: i) short-range con- 
tributions associated with momentum scales of order A^^^, encoded in the Lagrangian parameters and 
high-momentum part of pion loops; and ii) long-range contributions associated with scales of order m^r 
and below, which appear in the low-momentum part of loops and can be explicitly computed. The 
arbitrary separation between Lagrangian parameters and high-momentum part of pion loops is con- 
trolled by the regularization scheme, and after renormalization observables are regularization-scheme 
independent. 

In this context, the short-range EDMs di appear at LO for all sources. Long-range contributions 
resulting from pion loops introduce a dependence on g^^ and The specific chiral order at which 
a long-range contribution arises depends on the chiral properties of the underlying source, as detailed 
above. For the ^-term, CEDM, and left-right four-fermion operators, there are contributions propor- 
tional to g^"* at both LO and NLO, but at LO only to the isovector nucleon EDM. For the CEDM 
and left-right four-fermion operators there are additional NLO terms proportional to gl^'' as well. The 
dependence of loops on the regulator, choice of renormalization scheme, and renormalization scale can 
be absorbed into a renormalization of the d^. The momentum dependence, and in particular the con- 
tribution to the Schiff moment, is finite and set by 2mTj. For other sources, the relative suppression 
of pion-nucleon couplings means smaller loop contributions. The momentum dependence becomes a 
higher-order effect and the scale of its variation is determined by Mqcd rather than 2m^. 

For the nucleon EDM one obtains to NLO: 



and 



dn — dn 



di = di{fi) + 6di{fi) - 



egA 



egA 
(27r)2F^ 



3m^ 4(AmAr)5 



rriN 



In 



m 



mi 



4m 



N 



Ami 



mt 



rriN 



Tim^, 



Am 



N 



(3.107) 



(3.108) 



where Ami 



2 

m ± 



m 



^0 is the (mostly electromagnetic) pion mass splitting and 



Sdi{fj,) 



-(0) 

egAgir 



lE + In 



47r/i^ 



m 



N 



(3.109) 
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in terms of the dimension of spacetime d, the renormahzation scale /i, and the Euler-Mascheroni constant 
7e ~ 0.577. 

The /X dependence of 6di{n) can be absorbed in di{fi). One cannot separate parameter and loop 
contributions in a model independent way, and since the same combination of loops and parameters 
appears when the nucleon is inserted in a nucleus, it is simplest to redefine do and di to represent, 
respectively, the full isoscalar and isovector nucleon EDMs. However, one expects no cancellations be- 
tween loop contributions that are non-analytic in the quark masses and thus m^, and short-range pieces, 
which are analytic In this case the magnitude of the non-analytic contributions at a "reasonable" 
renormahzation scale serves as a lower bound for the redefined di. For n = itln, we expect for 6 [61] . 
qCEDM [To] and left-right four-fermion operators [22] , 

\do\ > 0.01 + 0.3 efm , ~ 0.1 [# + 0.03 #] efm . (3.110) 

Relying on the arguments leading to Eq. 
from which the current bound on 6 arises. For the other sources one currently has to rely on NDA or 
model-dependent estimates, as we discuss below. 



3.70) for for the 6'-term we obtain {d^l > 2 ■ 10 6 e fm. 



3.4 Lattice QCD 

Lattice QCD holds the promise of providing the various matrix elements for all relevant CPV mecha- 
nisms. However, to our knowledge, calculations have focused on the nucleon EDM for the ^-term, the 
isoscalar scalar couplings gs''^\ and the and tensor form factors. In the former case, the most recent 
published computations of the ajy date back nearly five years or more [751 ESI EH ESI ES EH] • Generally, 
these computations have followed one of two approaches: (a) computing the shift in the nucleon energy 
in the presence of an electric field, or (b) computing the nucleon electric dipole form factor by expanding 
to leading non-trivial order in 6. These calculations are carried out at unphysical values of the pion 



mass. The nucleon EDFF discussed in Section 3.3 provides in principle the tool to extrapolate results 



to smaller momentum and pion mass. Here, we summarize the most recently reported computations 
for each approach. 

The most recent computation of the first method has been reported in Ref. [79] . Using a two- flavor 
dynamical clover action, the authors considered the ratio of spin-up and spin-down nucleon propagators 

R,{E,t;e)=^-M^=[l + 0{e)]exp[-ar,eEt] , (3.111) 

where the subscript a on No- denotes spin, t denotes the time, and E gives the magnitude of the electric 
field along the z-direction. The electric field is introduced via a replacement of the gauge link variables 
Uk{x) in the Dirac- Wilson action as 

Uk{x) ^ e«'^'=*t/fc(a;) , (3.112) 



where Qq is the quark charge and k labels the direction. For this computation, the replacement (3.112) 



was applied only to the valence quarks; so-called "disconnected" insertions of the electric field on the 



sea quarks that enter through the quark determinant have not been included A corrected ratio R^°" 
was used to minimize the effect of insufficient statistics associated with vanishing E and/or 6. One then 
has 

'Rl°"{E,t-l;e 



2aN9E = In 



RT"{E,t- 



(3.113) 



^It is often conventional to retain only the non-analytic terms from loop computations and absorb all analytic terms 
into the parameters. 

^"in the limit of degenerate sea quarks in three-flavor QCD, the disconnected contribution is identically zero due to the 
vanishing trace over the quark charges |79) . 



24 



0.08 ■ 



0.04 - 



-0.04 - 



-0.08 - 



CP-PACS,N=2c 
CP-PACS,N=2c 
RBC,N=2DW,f 
QCDSF, N =2 do 
Current algebra 



0.2 



_1_ 



04 



_±_ 



0.6 

mAoeV') 



0.12 



0.08 



0.04 



0.2 



CP-PACS,N=2c 
CP-PACS,N=2c 
QCDSF, N =2 do- 



_L_ 



0.4 



_±_ 



0.6 

mAoeV') 



Figure 2: Lattice computation of nucleon EDMs induced by the QCD ^-term. The pion mass 
squared dependence of dn (left) and dp (right) obtained using various approaches. Square 
symbols denote the results in external electric field method in Nf = 2 clover fermion [79] , 
and circle symbols denote one in form factor method [HI] with same gauge configurations. 
Red bar denotes the bound of EDM in Nf = 2 domain- wall fermion in [77|, and diamond is 
a result from EDM form factor of imaginary 9 method quoted in [80]. Note that the error 
bar of diamond symbol may be an underestimate due to large systematic error associated 
chiral symmetry breaking of clover fermion. The triangle symbol is model estimate in current 
algebra. 



A 24^ X 48 lattice with /3 = 2.1 and lattice spacing a ^ 0.11 fm was employed, where the latter is set 
by the p-meson mass rrip = 768.4 GeV. Results for a„ obtained with a lightest quark mass corresponding 
to m,r = 0.53 GeV are shown in Fig. ^ using 9 = 0.025 and E = 0.004/a^. The corresponding values 



are quoted in Table 15 



The authors also studied the dependence of on the light quark mass to determine if this coefficient 
vanishes in the chiral limit as required. Results were obtained at = 1.13, 0.93, 0.76 and 0.53 GeV. 
Results for the neutron are indicated in Fig. [2j It is apparent that the computation does not exhibit the 
correct chiral behavior. The authors conclude that this situation is likely due to the explicit breaking 
of chiral symmetry by the Wilson-type quark action and the relatively large value of the lightest quark 
mass used. As the authors also emphasize, obtaining a significant, non-vanishing signal for the nucleon 
EDM does not appear to require the presence of appropriate chiral behavior. 

The most recent computation utilizing the form factor method has been reported in Ref. [HO]- The 
computation was performed by rotating 9 into the quark mass matrix and taking it to have an imaginary 
value: 

9 = -0 , (3.114) 

with 9^ being a real number. Simulations were performed using the Iwasaki gauge action and two-flavors 
of dynamical clover fermions with /3 = 2.1, a ^ 0.11 fm (again set by mp), rriT^/mp 0.8, and several 
values of the imaginary vacuum angle: 9^ = 0, 0.2, 0.4, 1.0, and 1.5. The EDM form factor was 
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Figure 3: Lattice computation of 6'-dependence of dn using the 
form factor method method [80j for 9^ = 0.2. Shown is the squared 
momentum transfer dependence at three mass parameters K = 
0.1382-0.1367 which correspond to ml = 0.3-0.85 GeV^. These 
are results in Nf = 2 clover fermion configurations. 



obtained from the ratio of three- and two-point correlators: 

R(t) = / 1 , (3-115) 

where t denotes the time co-ordinate for the insertion of the vector current J^, t' gives the time for the 
nucleon "sink", and p {ff) gives the nucleon momentum before (after) the vector curren insertion. 

Results at vanishing momentum transfer were obtained using two different extrapolation methods: 
(a) employing a dipole ansatz for the g^-dependence of the form factor and (b) assuming the EDM and 
Dirac form factors have the same g^-dependence and utilizing the latter (see Fig. |3]). Both methods 
give consistent values for the EDM. Taking 



d% = ^^ + c[er (3.116) 

and using the coefficient of the linear term to define the EDM, the authors obtain the results indicated 
in Table 15 The results agree with those of Ref. [79] (electric field method) within error bars. 

In addition to the direct computations of dN, lattice QCD results provide input for the determination 
of A(o) via Eq. ( |3.70 ) and for the g'g^'' ■ As discussed above, values of may be inferred from lattice 
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computations of a-i^N and (Am7v)g- Alternately, one may obtain g^^j, from direct computations of the 
charge changing scalar and tensor form factors[63] via isospin rotation. Taking into account the factor 
of two difference in normalization of these form factors, the preliminary lattice values quoted in Ref. [63] 
imply 

c/lj^^(MS,/i = 2GeV) = 0.4(2) (3.117) 
^^^^(MS,/i = 2GeV) = 0.53(18) . (3.118) 

The computation of g^p was obtained using two different gauge field ensembles with pion masses in 
the ranges 390 < < 780 MeV and 350 < < 700 MeV, respectively. A chiral extrapolation 
was performed assuming a linear dependence on m^. The value for g^^ was derived by combining 
RBC/UKQCD and LHPC results, with a chiral extrapolation based on HB^PT results. A comparison 
of the value for g^^"* with a result obtained using {AmN)q is given in Table 
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below. 



3.5 QCD Sum Rules 

In recent times, the most widely quoted hadronic computations of dn and the g^ rely on the method of 
QCD Sum Rules (QCDSR). (For an extensive review in the context of EDMs, see Ref. [H]; see Ref. |82] 
for a more recent discussion.) This approach entails computing hadronic correlators at large virtuality 
where the operator product expansion (OPE) can be rigorously applied and matching the result onto 
a phenomenological ansatz for the structure of the correlator at lower virtuality. The reliability of this 
matching is improved by performing a Borel transform to the OPE and phenomenological forms for the 
correlator. 

In the present instance, the relevant correlator n(Q^) involves two nucleon sources //at, 

Ii{Q^) = i j d^xe'^-^ (0|T{r/;v(a;)r/jv(0)}|0) , (3.119) 

where rjNix) contains combination(s) of quark field operators that carry the nucleon quantum numbers, 
where = —q"^, and where the Dirac indices on 11 have been suppressed for simplicity. In general, one 
uses a linear combination of two sources, 

VN = Vi + PV2 , (3.120) 

where 

?7i = 2tabcd^C^^Ubdc and = 2eabcd^Cuh^5Uc (3.121) 

with the subscripts a etc. denoting color. One computes n(Q^) in a background that contains photon 
and pion fields as well as the TVPV interactions introduced above and identifies various Lorentz struc- 
tures that are invariant under chiral rotations and that exhibit the appropriate spacetime symmetries 
associated with the EDM or ttNN interactions: 



EDM ^ 

vrA^A^ i , (3.122) 



is 



where F ■ a = F^^a^^. The corresponding phenomenological ansatz for the correlator at low virtuality 

n(Q2)Pheno ^ _ l/^(g2) |^ _ ^ \fM)i + " " " , (3-123) 



27 



with 



UQ^-) ^ (3.124) 



The dependence of the fk{Q^) on the quantities of interest here {dn, gi''^) arises from the first (double 
pole) term on the RHS of Eq. ( |3.124 ). As we discuss shortly, it carries an important universal dependence 



on the parameter A. The single pole term represents contributions associated with transitions between 
the neutron and excited states. Its strength, parameterized by the quantities Ad,.„ is generally unknown, 
leading to one source of theoretical uncertainty. The continuum terms Bd^n are generated by transitions 
between excited nucleon states and are also generally unknown. 

For the OPE evaluation, one performs all possible contractions of the q and q fields in the sources 
TjN and f/AT, leading to an expression in terms of the quark propagators evaluated in the presence of 
a photon and CPV background. The OPE gives the quark propagators in terms of various Wilson 
coefficients times condensates, viz. 

S{x) = + x(x)x(0) + ^(x)^P^°*°° + ^(x)isiuon ^ . . . ^ ^3^^25) 

where S{x)^^^ contains the free quark propagator, S{xY^'^°^'^^ and 5'(a;)^§^"°'^ contain dependences on 
the photon and gluon field strength tensors, respectively, and the x(^)x(0) carries dependences on the 
condensate (qq) as well as tensor condensates in the presence of the photon background that can be 
related to (qq) through various susceptibilities x, i^, and ^: 

{qcyi,uq)F = QqX F^,u{qq) , (3.126) 
gs{qGf,uq)F = Qqfi Ffiuiqq) , (3.127) 
2^s(gG'^.75g)F = iQq^ F^,{qq) . (3.128) 

After including the CPV interactions in the background, one obtains the dependence of x{^)x{^) on {qq), 
the susceptibilities x, f^, and ^, and the quantities 9, dq and dq. A detailed expression for the resulting 
correlator n((5^)°™ goes beyond the scope of this review but can be found in, e.g., Refs. j83 | IM j l85 | 182] . 

Applying the Borel transformation to both T{{Q'^Y^'^^" and n(Q^)^^^ and matching the coefficients 
of the relevant Lorentz structures then allows one to obtain dn and the g'i^'^ in terms of the CPV 
parameters, susceptibilities, {qq) condensate, Borel mass M, and phenomenological parameters A, Ad, 
etc. For example, for the neutron EDM one has [1] 



M 



4 



X'dnrriN - AM' = -0 {qq) ^e"^"/^' , (3.129) 
with 

e = {AQdnidPd - QurriuPn) xO + (4rfd - d^) + {k - ^/2) (iQddd - Qudu) , (3.130) 

and with the pq carrying a dependence on ratios of the quark masses. In the absence of a PQ mechanism, 
one has 6 = 6. As discussed above, in the presence of the PQ mechanism the other QCD CPV 
interactions, such as the CEDM, lead to a shift in the vacuum angle due to their effect on the axion 
potential. In this case, one must take 6 = 6ind, with 6i^d being the shift due to the additional axion 
potential contributions. 

The most recent results [HI [82] for the dependence of dn and the gi!'^ are indicated in Tables [l5 



19 We concentrate first on the dependence of dn on 6. Notably, the magnitude of obtained from 
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the QCDSR computations are one to two orders of magnitude smaller than those obtained using lattice 
calculations. Moreover, the most recent QCDSR determination of a„ [S2] is nearly a factor of six smaller 
in magnitude that the earlier work of Ref. [11]. This difference results, in part, from a different value 
of A used in the two analysis. The authors of Ref. [H] utilized the smaller value for A obtained directly 
from QCDSR studies, while Ref. [82] employed a value obtained from a lattice computation of proton 
decay matrix elements and isospin symmetry, which is roughly two times larger. 

The foregoing technique also provides the dependence of dn on the quark EDM and chromo-EDMs. 
Again, the analyses in the two most recent computations [HI 182] are similar, differing primarily in their 



extraction of A. Illustrative results are given in Tables 16 and 17 for the case when one assumes the 



presence of PQ symmetry. To our knowledge, the only QCDSR computations of the dependence of the 
g^^ on 6, dq, Cq and four-quark operators have been reported in Refs. [861 HI]. As discussed above, one 
may extract the ^-dependence of g^^ using chiral methods, so we do not quote an additional QCDSR 
result for this case. Moreover, the contribution of 9 to g^^ requires an additional power of ruq (or m^) 



as indicated in Eq. (3.73), so its impact will generally be negligible except for systems in which the 



effect of (7^°^ vanishes. Consequently, we do not list any entry for A(i). 

For the contributions of other CPV operators, the situation is more subtle. In the case of CEDM 
contributions to the g^\ the constraints of chiral symmetry imply the presence of two canceling contri- 
butions: 

{N\ {dq ggsCr^.G^^'q - mlqq) \N) (3.131) 

where the first term arises from the PCAC commutator term and the second is generated by the pion 
pole contribution. This cancellation renders the computation susceptible to theoretical uncertainties, 
particularly associated with the choice of (3 in the nucleon sources. For example, taking /3 = 1 yields 
a vanishing contribution through NLO under the assumption of pure valence quark dominance. On 
the other hand, choosing /3 7^ 1 yields a non- vanishing result. Assuming that the double-pole term in 



Eq. (3.124) dominates, one obtains 



^i"-^3^-'(«(^»--^)' 

where at leading order the -Ffc(/3 = 1) = and -Ffc(/3 = 0) = 1. For the latter choice, one obtains a 
generically stronger CEDM sensitivity of glr'^ compared to that of glr'K Going beyond LO, including 
uncertainties associated with additional condensates that consequently appear, the choice of Borel mass 
M and variations with /3, the analysis of Ref. [86] yields a "best value" and range for gl^^ as well a 
broad range but no best value for ^1°^ quoted in Table 17 

Providing robust computations of the three-gluon and four-quark operators to the neutron EDM 
and the cjn^ is even more challenging, as one encounters additional unknown condensates as well as the 
presence of infrared divergences at lower order in the OPE than for the other sources of CPV discussed 
thus far. It is possible, however, to estimate the contributions to dn by relating the EDM to the nucleon 
magnetic moment through a CPV rotation of the nucleon wavefunction. One then has, for example 

9gsml ( M'^\ /v\2 



where fXn is the neutron magnetic moment and where the other pref actors in Eq. (3.134) arise from a 
QCDSR evaluation of the correlator of two nucleon currents in the presence of the CPV three-gluon 
interaction. The latter evaluation is used to determine the CPV rotation needed for rela ting dn and fir, 
Taking M/Ajr = 2 and gs = 2.1 yields the estimate of the coefficient (3^ given in Table 
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3.6 Quark Models 

The constituent quark model (CQM) has proven remarkably successful in accounting for a number 
of static properties of the lowest-lying baryons, most notably their magnetic moments. In the latter 
instance, one assumes each constituent quark posses a distinct magnetic moment that is proportional 
to its spin, 

/iQ = 2fiQSQ , (3.135) 

and computes the nucleon magnetic moment using the appropriate spin-fiavor-color-spatial nucleon 
wavefunction, resulting in 

A^P = ^ - f^o] and = ^ ['^'"^ ~ ■ (3.136) 

To the extent that the constituent up- and down-quark magnetic moments differ only by the overall 
charge of the quark, one obtains for the ratio of nucleon magnetic moments fip/fin — —3/2, in close 
agreement with the experimental value. 

In the case of the EDM, the dimensionless Wilson coefficients Cg^ and 6q (or the equivalent di- 
mensionful quark EDMs dq) correspond to the EDMs of the current quarks of QCD rather than the 
constituent quarks of the quark model. Nevertheless, one may make the bold ansafe that the coefficients 
of the constituent quark operators are given by the corresponding coefficients for those of the current 
quarks after appropriate RG running from the scale A to the hadronic scale that introduces the Kg 
factor in Table [5l 

dg = 2dQSQ with dg = dgiA^) = Kgdg{A) + " " " , (3.137) 

where the last relation may just as easily be expressed in terms of Cq^ or 6q and where the indicate 
contributions associated with operator mixing. In this case, the computations of the nucleon EDMs 
proceed as in the case of the magnetic moments, leading to 



dp = - [Mu — do] and dn = - [4c//) — du] ■ (3.138) 



Thus, one obtains 



pt = p; = -^Pl = -4pJ = - . (3.139) 



4 
3 

An alternate approach, first proposed in Ref. [HE] , is to retain the identity of the quarks as partonic 
degrees of freedom and relate the nucleon matrix elements of qa^^^qF^^^ to those of the quark axial vector 
currents q'^^.'^^q that contribute to the nucleon spin: 

l{N\qi,l5q\N) = (Aq)^ S, , (3.140) 

\{N\qa^.^5q\N) = {Aq)^Na,,N, (3.141) 

so that the nucleon EDM is given by 

dN= J2 ^^(^^)7V ' (3-142) 

q=u,d,s 

leading to a "parton quark model" (PQM) prediction 

p^ = (Ag)^ . (3.143) 

Information on the (Ag)^ can be obtained from a number of sources, including the determination of 
the axial vector coupling that enters neutron decay, 

gA = (Au)^ - (Aw)„ , (3.144) 
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and studies of both inclusive and semi-inclusive polarized, deep inelastic leptoproduction, 

(An)^ = (Arf)„ = 0.746 , (Arf)^ = (Au)^ = -0.508 , (A.)^ = (A.)„ = -0.226 . (3.145) 

The CQM and PQM fail to provide guidance for the quark CEDM or three-gluon operator contri- 
butions to the nucleon EDMs. They similarly offer no reliable means for estimating the contributions 
of these operators to the g^^ . In these cases, one might complement quark model estimates with the 
NDA discussed in Section l3^ 



3.7 Saturation Methods 

Over the years, estimates of four-quark matrix elements have often been obtained by assuming certain 
intermediate states dominate or "saturate" the dynamics. Here, we illustrate the application of the 
saturation approximation to estimate two classes of matrix elements: (a) contributions to glr^ generated 
by the operator Q^ud and (b) the values of the TVPV fou r-nuc leon operator coefficients C'i_2. 



Starting with gl^\ we recall that the first term in Eq. (|3.89|) generates a contribution to g^''' ■ When 



the scale 1/A^ is included, one expects ^l^"* ~ ImC^ud-^x-^-^/-^'^ noted earlier. An explicit evalu- 
ation can be made using factorization and vacuum saturation, partial conservation of the axial cur- 
rent (PCAC), and the vrA^ cr-term. Vacuum saturation in this context amounts to first relating the 
(A^7r2| 5*3 ® 5*4 |A^) to the crossed matrix element {n2 \ S3 ® 5*4 \NN), inserting a complete set of states 
between the S3 and 5*4 bilinears, assuming the dominant contribution arises from the vacuum, and then 
uncrossing the N. One then obtains 

{Nn2\ S3 ® ^4 |A^) ~ (A^l ^4 |A^) (vTsI ^3 |0) . (3.146) 
Now, the matrix element of 6*4 is related to the pion-nucleon a-term as 

(N\S^\N) = {N\uu + dd\N) = — , (3.147) 

m 

while the second matrix element can be evaluated by taking the divergence of 

Mp)\q'^Yl5q |0) = -iF^p^ (3.148) 



and taking the pion on-shell, leading to 



(vr3|53|0) = ^^ . (3.149) 
m 



Using Eqs. (3.146 3.149) and including the coefficient of the four-quark operators appearing in Eq. (2.29) 
leads to 

(1) /T ^ \ ( F-k<^-kN 



t^^ = {\^C,^d)[^-j^) y-^) , (3.150) 
which is of order k^F^,! as advertised. Using a-„N ~ 45 ± 6 MeV[53l [M], m ^ 3.85 MeV [21] gives 

# = (3.3 X 10-5) X (ImC^,,) (^)' . (3.151) 

The other approach is to assume the saturating states are the lowest available single-meson states. 
For example, a meson of mass mm 3> m^r that can be exchanged between two nucleons gives rise to a 
potential of range ~ 1/mm, which is short compared to typical nuclear distances ~ 1/m^ ^ Xjmra- In 
an expansion in powers of m.„/mra, such a potential can be replaced by contact interactions with an 
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increasing number of derivatives. In first order in the relevant TVPV parameter, the meson couples 
through a TVPV coupling to one nucleon and a P-, T-even coupling to the other, resulting in TVPV 



contact interactions such as the Ci,2 terms in Eq. (3.38). 



In this context, the mesons that have been considered are the lightest: the pseudoscalar i] and the 
vector mesons V = p,LJ [891 EOl [91] , with TVPV interactions given by 

ClY"" = VN + 4^V3) N + ^[N {-g(^^ + -gWr,) ^a^^^,^^N + H.c] 

+ { [^Wr ■ + -gf^pl + (3r3p° - r • pj] la^'^i.d^N + H.c.} . (3.152) 



N 

They lead to 

1 / -(0) -(o)\ 
^ rriN y ml j ' 

Ci'' = -^(^f , (3.154) 

where and my are meson masses, and Qr^NN and Qvnn are P-, T-even meson-nucleon couplings, 
respectively the axial coupling of the eta and the vector coupling of the vector meson. This type of 
meson exchange also produces other contact interactions [91j, which are, however, expected to be of 
higher order for CPV sources of dimension up to six. 

Since my ~ and there is no reason for (jy^ /glr'^ to be particularly big or small, the size of vector- 
meson contributions is comparable to the NDA expectations, with some suppression coming from the 
numerical smallness of the P-, T-even rho-nucleon vector coupling QpNN — 3.2 [021113] compared to the 
analogous pion-nucleon coupling 2mNgA/ F.„ ~ 13.5. (In contrast the same ratio for the omega is close 
to 1 [921 |93].) For the eta meson, the enhancement due to the relatively light mass is offset by the 
relative smallness of QrjNN — 2.24 [M]. Obviously the limitation of this type of saturation is that there 
are no firmer estimates of the TVPV couplings and cjy^ than for the 6*1^2 • 



3.8 Hadronic Matrix Elements: Discussion 

While there exists a solid body of work devoted to matching the ^-term and dimension six operators 
onto hadronic quantities, there clearly exists considerable need for further advances. In what follows, 
we attempt to provide a sense of the present range of theoretical uncertainty in sensitivity of various 
hadronic quantities to the underlying operator coefficients. To that end, for each sensitivity coefficient 
(«„, A(o), Cn; ^^c.) we provide a "best value" and "reasonable range". The importance of attempting to 
quantify the theoretical uncertainty is two-fold. First, when using EDM search limits to derive bounds 
on the underlying parameters such as 6 or 6q, previous studies have often included only the experimental 
uncertainty while relying on a single hadronic computational framework. As a result, the quoted bounds 
may be unrealistically stringent. Second, we anticipate that hadronic structure theorists involved in 
the relevant computations will find our benchmarks helpful in setting goals for future refinements. In 
short, it is useful to know where and to what level reductions in theoretical uncertainties are called for. 

Before proceeding, we make an important caveat. In comparison with the analysis of experimental 
uncertainties, the task of assigning theoretical error bars entails a greater degree of subjectivity. Con- 
sequently, the benchmarks we provide below should be taken somewhat impressionistically rather than 
as quantitatively robust. Nevertheless, we believe they offer reasonable guidance as to the present level 
of uncertainty as well as quantitative targets for further refinements. 

With these considerations in mind, we now discuss specifics. In setting the best values and ranges, 
we will use the considerations based on chiral symmetry/NDA as well as the spread of current theoretical 
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computations as guides. Chiral symmetry and NDA is particularly helpful in determining if the results 
of specific computation are anomalously large. While one expects the chiral arguments to be realistic 
up to a factor of a few, a computation that yields a result an order of magnitude larger would likely 
be open to question. On the other hand, it is quite possible that specific dynamics can suppress a 
given quantity by more than a factor of a few. In general, our best values will be somewhere close 
to the mid-point of the range of explicit computations and close to the magnitude indicated by chiral 
symmetry/NDA. When assessing the range of explicit computations in a given framework, we will be 
rather inclusive, except when a given computation seems to be particularly anomalous with respect to 
the chiral/NDA expectations. Hence, our reasonable ranges will be roughly consistent with the spread 
of explicit computations and the factors of a few variation one might expect with respect to chiral/NDA 
arguments. 

Table |6| then, gives these benchmarks. Generally speaking, we see that the quantities with the most 
narrow ranges are: 

(a) a„, the sensitivity of dn to 6: QCD sum rule computations are quite in line with chiral symme- 
try/NDA expectations. Moreover, a naive scaling of the lattice results with would imply a 
value close to our best valu^^ Note that since the ^-term does not run, there exists no uncertainty 
associated with an incomplete analysis of RG evolution, in contrast to several of the dimension 
six operators. 

(b) A(o), the sensitivity of g^^ to 6: the use of chiral symmetry and lattice results for {AmN)q provide 
a relatively model-independent result. Reduction in the errors on (AmAr)^ and the light quark 
masses from the lattice will lead to a corresponding narrowing of the theoretical range on this 
quantity. 

(c) g^g \ isoscalar scalar form factor that governs in part the sensitivity of atomic and molecular EDMs 
to the combination of coefficients lm{Ciedq — C^^lu)'- a model independent value is obtained from 
a-n-N and the average light quark mass, rh. To the extent that the lattice uncertainties on these 
quantities are robust, one has a relatively narrow range for the isocscalar scalar form factor. 

(d) Pn, Cn)y the sensitivity of dn to the quark EDMs: results of explicit computations do not vary 
considerably from expectations based on either chiral symmetry/NDA or the quark model. It is 
worth emphasizing, however, that a complete analysis of the RG evolution of the quark EDMs 
from the weak to hadronic scales, taking into account mixing with the CEDM and four-quark 
operators, has generally not been carried out. 

Although further reductions in the uncertainties associated with these quantities would be welcome, 
we do not consider them to have the greatest urgency. Those seemingly most theoretically fraught are 
the sensitivities to the CEDM, three-gluon operator, and four-quark operators. 

(e) (/S^*^, Pn, Cn)? the sensitivity of dn to the quark CEDMs: Here we take as best values the average 
of the existing QCD sum rule results, which in the case of the d-quark CEDM is equal to the 
chiral/NDA expectation. The ranges here are rather broad, spanning an order of magnitude. 
Moreover, as in the case of the quark EDMs, a complete implementation of the RG evolution 
that includes mixing with the four quark operators remains to be performed. We also note 
that recent studies of the CEDM contribution to the p-meson EDM using the Dyson-Schwinger 
approach [95] raise further questions about the reliability of matching of CEDMs onto hadronic 
quantities. Although the p-meson EDM is not of experimental interest, the relatively simplicity 

^^Onc should note, however, that the lattice computations to date do not necessarily manifest the expected chiral 
scaling in other related observables, so the result of naive scaling may be a coincidence. 
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of the p-meson bound state makes it a useful "laboratory" for various hadronic matrix element 
computational methods. The results obtained with the Dyson- Schwinger framework imply that 
the p-meson EDM is an order of magnitude less sensitive to the CEDMs than one would infer from 
the corresponding QCD sum rule computation [96]. In contrast, both approaches yield similar 
sensitivities to the quark EDMs. Should a similar situation emerge for dn-, one would need to 
further inflate the theoretical uncertainty on the (/S^*^, p^, Cn)- 

(f) {ifi)-, V{i))> the sensitivity of the glr^ to the quark CEDMs: Here the situation is even more 
uncertain. To our knowledge, only two computations have appeared to date. The ranges quoted 
in Ref. |86] is consistent with the magnitude expected from chiral symmetry/NDA, but in the 
case of the contribution to glr'^ incompasses zero as well. In the case of the g^^ sensitivity, the 
range in Ref. |86] is narrower, and we have no present rationale to expand it, but the dearth of 
analyses and the range for the glr^^ sensitivity should give one pause. 

(g) and 7^^^, sensitivity of dn and the g^^ to the three gluon operator: For this case, we posses 

a dearth of information. The central value for given in the QCD sum rule work of Ref. [87] 
that is often quoted elsewhere is an order of magnitude smaller than the chiral symmetry/NDA 
expectations. Consequently, we take a broad range for this parameter. For the 7^*^^, we have only 
the chiral/NDA expectations, and, thus, employ the "factor of a few" criterion for this range. 
Note that the three-gluon operator is multiplicatively renormalized, so the theoretical uncertainty 
is associated entirely with the hadronic matching computations. 

(h) Four quark operators: hadronic matrix element computations for these operators is, if anything, 
even less advanced than for the three-gluon operator. Apart from issue of RG evolution, explicit 
computations are few and far between. The factorization computation used to match Q^ud onto 
(jn'^ gives a smaller value for 7^'^"'^ than the chiral/NDA expectation. Computations for the other 

operators Q^J^^ have been carried out using a combination of the quark model, factorization, and a 
relativistic meson loop approach[97]. However, the meson loop computation utilized in that work 
is not consistent with the EFT power counting embodied in HB^PT and, thus, may overestimate 
the magnitude of the matrix element by an order of magnitude. Consequently, we are reluctant 
to use the results in that work for guidance. Instead, we start with the chiral/NDA estimates 
and give an order of magnitude spread based on the present dearth of consistent computations, 
implementation of RG running, and comparison with the factorization estimate. 

We do not include in Table [g] the semileptonic form factors apart from g^g'^^- The manifestation of 
the pseudoscalar and tensor form factors in atoms and molecules is suppressed by several factors. The 
pseudoscalar interactions are higher order in the HB^PT expansion, while the tensor charge does not 
receive a nuclear coherent enhancement. We note, however, that the EDM of the diamagnetic atom 
i99j^g ]-^g^g roughly an order or magnitude greater sensitivity to the tensor eg interaction than it does to 
the scalar interaction. The value of gf^'^^ in this case is, thus, considerably more significant than in the 
paramagnetic systems. 

Two implications should be drawn from our theory uncertainty estimates. First, for BSM scenarios 
in which the CEDM, three-gluon, and/or four-quark operators have significant Wilson coefficients, it 
will be particularly important for any global analysis to include the rather sizeable uncertainties at 
the hadronic level. Second, a concerted effort to refine the hadronic computations and reduce the 
uncertainties is clearly in order. We hope that our delineation of these best values and ranges will spur 
future efforts in this direction. 
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4 Nuclear, Atomic, and Molecular Scales 



Composite systems are often the most amenable to experiments. But the EDMs of composite systems 
reflect few- or many-body dynamics as well as the fundamental source of CP violation and QCD. For 
some experiments, completely ionized light nuclei are useful, and for these cases we must employ few- 
body methods to relate the nuclear EDMs to, e.g., 9. In heavier neutral systems a new phenomenon 
is important: the shielding of the EDMs of constituents of one charge {e.g. protons in the nucleus) by 
those of the other (electrons). The transmission of CP violation through a nucleus into an atom must 
overcome this shielding and its effectiveness in doing so is expressed by a nuclear Schiff moment, which 
we define shortly. 

We begin by considering EDMs of light systems, potentially important for storage-ring experiments, 
and then move to heavier systems, useful in experiments on immobilized atoms or molecules. 



4.1 Light Nuclei 

In addition to the continuous improvement in experiments on neutral systems, a new, exciting prospect 
is the direct measurement of the EDMs of charged particles in storage rings [981 EH] • When a particle 
moves in an electric and/or magnetic field, its spin will precess at a rate that depends not only on the 
magnetic dipole moment but also on the EDM. The best bound on the muon EDM |lUUj comes, in fact, 
as a by-product on the BNL g — 2 experiment. It can be expected that dedicated experiments in rings 
with optimized parameters will allow sensitive probing of the EDMs of light nuclei. For example, it has 
been proposed [98l 199] that for dp a sensitivity of 10~^^e fm can be achieved. Similar sensitivity could 
be attained also for the deuteron (^H nucleus) and helion (^He nucleus) EDMs, dd and dh respectively. 
The triton (^H nucleus) EDM, dt, might be accessible as well. 

From a theoretical perspective, the EDMs of light nuclei can be calculated with relatively small un- 
certainty originating in the P-, T-even strong interactions, as essentially exact calculations are possible. 
Moreover, with an effective field theory approach based on HBxPT we can treat the nucleon and nuclear 
EDMs on the same footing, and explore the sensitivity of nuclear EDMs to different combinations of 
TVPV hadronic interactions than that which appears in the nucleon EDM. In particular, it has been 
argued that the deuteron EDM has some sensitivity to the CP V source |101t 1102] , and that a combined 
measurement of d^, dp, dd, d^ and dt could be used to disentangle the various sources [103j. The reason 
for this is the different relative strengths of the various couplings at LO in HB^PT jlO^ I1U3] . which 
are rooted in the different chiral symmetry properties of the various sources, as discussed in Section 



3.2[ Likewise, experimental access to other TVPV moments, such as the deuteron magnetic quadrupole 
moment, would be very useful as well for separating sources |102[ 1104] . but it does not look feasible in 
the near future. 

In a nucleus with A nucleons, certain P-, T-even inter-nucleon interactions need to be resummed 
in order to produce a bound state and its associated wavefunction \^ a)- (For a review of nuclear 
EFT, see Ref. |1U5] .) A nuclear EDM dA arises from the average with such a wavefunction of two 
TVPV mechanisms: i) the TVPV electromagnetic current Jtvpvi whose one-nucleon component is the 
nucleon EDM; and ii) a. combination of the TVPV potential Vtvpv and the P-, T-even electromagnetic 
current J^rp. We follow here Ref. |103] . which we refer the reader to for more details. Because TVPV 
interactions are so tiny, we can write in first-order perturbation theory 



dA = {^ 



D 



TVPV 



^4 ) + 



D 



PT 



^a) + c.c.) . (4.155) 



The electric operators are obtained from the corresponding time-component electromagnetic current 
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through D = ilim^^Q'VcjJ^{q). At the one-nucleon level 



D 



(1) 

PT 



i=l 



Si) 



i=l 



(4.156) 



in intrinsic coordinates, which obey ^ 



A ^ 



0. The more complicated two- and more-nucleon currents 



are expected to be generically less important, although this is not always true, as discussed below. The 



first term in Eq. (4.155) represents the contribution of the individual nucleons to the nuclear EDM, 



as well as the contribution from TVPV many-body currents. The second term in Eq. (4.155) is the 



contribution of the parity-admixed wavefunction a)^ obtained from the TVPV potential via 

{E - Hpt)\^ a) = Vtvpv\^ a) where {E - H pt)\'^ a) = ^ , 



(4.157) 



with HpT being the P-, T-even Hamiltonian. 

Because of its long range, one-pion exchange (OPE) has long been recognized as a potentially 
important component of the TVPV two-nucleon (NN) potential |106[ 11071 1108] , and expressed in terms 
of the three non-derivative pion-nucleon couplings in Eq. (3.38) [38] • In the literature, this potential 



is sometimes supplemented by the single exchange of heavier mesons, with the t] |109] , rho |110] , and 
u) |110] being most popular. Allowing sufficiently many couplings of these mesons to nucleons one can 
produce |lllj the most general short-range TVPV A^A^ interaction with one derivative |112j . Although 
a derivative expansion is justified on general grounds, the relative importance of terms with various 
ranges, spin/isospin structures and number of nucleons depends on the TVPV source. The potential 
obtained from HB^PT, which accommodates the most important effect of each source, is discussed in 
detail in Ref. [STj. In configuration space, the A^A^ potential corresponding to the Lagrangian (3.38) 
reads 



TVPV 



' — a 



U) 



) + 



• vr(r,w 



(4.158) 



where rij = rj — r*,- is the relative position of the two interacting nucleons and Y{r) = exp(— m7rr)/47rr 
is the usual Yukawa function, so that 



Vy(r,: 



47rr| 



1 + m^rij) exp{-m^rij) 



(4.159) 



The short-range interactions can be thought of as accounting for heavier-meson exchange, as discussed 
in Section |3.7[ For all sources, few-body potentials are expected to generate smaller contributions. 



except for the left-right four-quark operator (2.29), for which the effects of a three-nucleon potential 
originating in a LO three-pion vertex [103j proportional to gli^^ remain to be studied. 

Equation (4.155) has been evaluated in the literature for A = 2,3, and the explosive growth in ab 



initio methods affords ways to calculate the EDMs of larger nuclei if needed. Most existing work employs 
TVPV one-meson-exchange potentials, with older references using simple P-, T-even wavefunctions and 
single-nucleon currents, and more recent ones, highly developed phenomeno logical P-, T-even potentials 
and even meson-exchange currents. In HB^PT the non-analytic behavior of the nucleon EDM in m^r 
and the dominance of OPE in nuclear observables can be accounted for simultaneously, with chiral 
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symmetry playing a central role. In principle full consistency can be achieved, but so far calculations 
are limited to phenomenological P-, T-even interactions for Hpx- At this stage leading contributions, 
with an uncertainty of roughly ~ mT^/MqcD ~ 20%, are sufficient, and for the most part one can restrict 
oneself to the one-body currents and two-body potential described above. For the ^-term the situation 
is more complicated because in nuclei with equal numbers of protons and neutrons, N = Z, the isoscalar 
component of the P-, T-odd potential (4.158) gives a vanishing contribution in combination with Dpi' 
(4.156) [106j. The latter is an isovector with a conserved third component, and can only contribute to 
the EDM if there is a parity-admixed component of the wavefunction that differs from the ground-state 
wavefunction by one unit of isospin. For dimension-six sources, where the leading A^A^ potential is not 
expected to be dominantly isoscalar, this is not of particular consequence. But, for the ^-term, g^'^ 
is the formally leading part of the potential, and in A^ = Z nuclei a non-zero result comes only from 
subleading parts of the A^A^ potential as well as two-body currents. 

Not all the seven parameters shown explicitly in Eq. (3.38) are important for every CPV source. 

_(2) 

shows, gn is expected to be small for all sources, and as a 



In fact, as the discussion in Section 



3.2 



,(0,1) 



consequence the EDMs of light nuclei should be described at LO in terms of the szx parameters (io,i, g^ 
and Ci^2 |103] . Results for the EDMs of light nuclei in terms of these six LO parameters are reviewed 
below and summarized in Table [7j The first two rows are a reminder that we have absorbed the loop 



contributions to the neutron and proton EDMs in cIq i, as discussed in Section 3.3 The potential-model 
dependence in the subsequent rows is not larger than ~ 25%, which is comparable with the LO HBxPT 
error. Exceptions are the short-range contributions from the Ci, which can only be considered order- 
of-magnitude estimates. As discussed below, for the tri-nucleon system there are disagreements in the 
OPE estimates of about ±50% in the values quoted. After discussing specific results for the deuteron, 
helion and triton, we cast them in terms of the ^-term and dimension-six sources. 



4.1.1 Deuteron 

The deuteron EDM has been investigated in the meson-exchange picture |114l I115[ 11161 IllH 1117] , with 
various degrees of sophistication in the treatments of the P-, T-even interaction HpT- Refs. |11H I117j 
found that differences in the EDM generated by TVPV OPE are rather small among modern high- 
quality phenomenological potentials. Ref. |103] uses the calculation scheme of Ref. |lllj to obtain 



wavefunctions |\E'2) and |\E'2) for the potential (4.158), in conjunction with HBxPT currents, for all 



CPV sources of dimension up to six. Ref. [50j also uses HB^PT and phenomenological potentials, 
but with slight different dimensional estimates for the various HB^PT contributions, in the particular 
case of the 6 term. All these references ignore relativistic corrections, which are absent from the 
phenomenological potentials they use. A fully consistent HB^PT calculation exists |102] . in which 
pions are treated perturbatively — a good approximation scheme in the loosely bound deuteron. In 
this case the LO P-, T- even potential is just a delta function, so this is an effective field theory extension 
of Ref. |116] . Results are consistent with Refs. |103l 150] . suggesting that a fully consistent calculation 
with non-perturbative pions will not deviate significantly from the results obtained so far. 

The simplest contribution to the deuteron EDM originates in the EDMs of its constituents. Since 
the deuteron has spin S = 1 and isospin / = 0, the nuclear matrix element of D!^ypy in Eq. (4.156) 



gives simply 2do. Since there is no reason to expect a cancellation with other contributions, 2do serves 
as a lower-bound estimate for the deuteron EDM. For the 6'-term in particular, using the long-range 



NLO contributions to the isoscalar nucleon EDM, Eq. (3.110), we expect [102j 



Irfrfl >3- 10-^^efm . (4.160) 

In agreement with the more general argument for N = Z nuclei, the one-body operator Dp], cannot 
bring the deuteron wavefunction back to a (mostly) ^Si wave once the isoscalar TVPV potential takes 
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it to ^Pi. In order for Dp], to yield a non-zero contribution, the parity-admixed component of the 
wavefunction has to be in the ^Pi state, to which only glr^ contributes. The corresponding nuclear 
matrix element has been calculated several times in the literature. As summarized in Ref. [50], there 
is agreement to better than 10% among modern NN potentials, and within 30% between them and a 
simple delta-function potential. When OPE through the isotensor pion-nucleon coupling g^^ is included, 
it gives small contributions , even when estimates about the small magnitude of g^r are disregarded. 

The expectation from HB^PT that two- and more-body currents give small contributions is corrob- 
orated by a model calculation . They can be neglected for all CP V sources except the 6'-term. For 
the latter, because the formally LO contribution vanishes, one has to go to NNLO. This brings in the 
same dependence on g^^ as for other sources. Additionally, a dependence on (7!°^ emerges through the 
sub leading potential and two-body currents, together with two P-, T-even isospin-breaking parameters: 
the quark-mass component of nucleon mass difference, {AmN)q in Eq- (3.64), which can be estimated 



from lattice QCD as we have done in Section 3.2 and the isospin-breaking pion-nucleon coupling 



{(3l/F^){^^7f3)NS^'N, with (3^ = (^(em^/M^^^), for which only the bound /3i = (0 ± 9) ■ lO'^ 011 [TT3] 
is known. 

These results are summarized in the third row of Table |7l 



4.1.2 Helion and triton 

There have been fewer calculations of the trinucleon EDMs. A pioneering calculation |118j of the 
helion EDM for the glr^ OPE (dominant for the ^-term) used an old phenomeno logical P-, T-even 
potential solved in the adiabatic approximation of the hyperspherical-harmonics method, and found no 
nuclear enhancement compared to the neutron EDM. The era of modern calculations began with Ref. 
|119] . when the nuclear wavefunction was calculated using high-quality P-, T-even potentials including 
the Coulomb interaction. A solution is found with the no-core shell model (NCSM) method, which 
employs a model space made from Nmax properly antisymmetrized harmonic-oscillator wavefunctions 
of frequency Q. At large enough N^ax, results for the helion EDM, which are somewhat larger than Ref. 
[IIS J where they can be compared, become independent of Q. Meson-exchange currents were neglected, 
as suggested by their smallness in the deuteron. For mesonic couplings of equal magnitude, OPE is 
found to be dominant over shorter-range interactions. Ref. |103] adapted this calculation to the TVPV 
ingredients from HBxPT, and calculated the EDM of triton for the first time. The two short-range 
interactions from Ci^2, which can be thought of as originating from, respectively, omega and rho meson 
exchanges considered in Ref. [119j , were regulated with Yukawa functions. Ref. \120\ used similar input 
but solved Faddeev equations instead. 

Calculations with various realistic potentials agree within ~ 25% for the nucleon EDM contributions. 
They give nuclear matrix elements of roughly equal magnitude for c/q and di, so that, as one might have 
expected, the helion (triton) EDM is mostly sensitive to the neutron (proton) EDM, 

For the contribution from the TVPV potential, results for triton are very similar in magnitude to 
those for helion. Both Refs. |119] and |120j find a spread of ~ 25% between different potentials, but 
they disagree by an overall factor of about two in isoscalar and isovector TVPV terms (and a factor 
five in the subleading isotensor component). The reason for this discrepancy is unclear at present, and 
it is a priority to resolve it. Additionally, in Ref. |103j the short-range contributions from the Ci were 
found to vary considerably with the explicit regulator mass from potential to potential (a factor ~ 5 in 
the cases studied). There is thus a much stronger potential dependence, and more solid numbers have 
to wait for a fully consistent calculation. 

These results are summarized in the fourth and fifth rows of Table [7| where for OPE we took values 
in between those of Refs. |119[ 1103] and |120] . Note that there is yet no estimate of the effects of 
a three-nucleon TVPV potential, which could be significant in the case of the left-right four-quark 
operator [22] . 
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4.1.3 Light nuclear EDMs: summary 



Using tlie results of Section [3} we can obtain the sensitivity of each EDM to the underlying CPV sources. 
The orders of magnitudes expected for these EDMs are given in Table |8] for each source. Recall that the 
Weinberg three-gluon operator (2.23) and the two four-quark operators (2.25) present already at the 
electroweak scale are chiral invariant. As a consequence they produce the same hierarchy of hadronic 
interactions and cannot be separated at low energies. For simplicity we use the shorthand notation 
{Ca,lmd^l;^l} ^ ImCk in Table [s 

The texture of this table underlines the argument |10H 11021 1103j that light nuclei act as a "chiral 
filter" for the various CPV sources. Of course, a measurement of the neutron EDM alone could be 
due to a 6'-term of just the right magnitude, or to any of the dimension-six sources, although if all 
dimensionless factors were equal, chiral-invariant sources or the "left-right" four-quark operator (2.28) 
would be favored because they require no chirality flip. Nuclear effects are most significant for the 6- 
term, CEDM and left-right four-quark operator. Just on the basis of orders of magnitude, we see that a 
large tri-nucleon EDM compared to a nucleon EDM would point to them as possibly dominant sources, 
while a large \dd\ compared to {d^l would be suggestive of just the CEDM and left-right four-quark 
operator. Bounds on light EDMs would provide tighter bounds on this physics than comparable bounds 
on nucleon EDMs. 

We can infer more information about CPV sources from Table [7] when we take into account that the 
relative importance of various pion-nucleon and short-range interactions is not the same for all sources. 
For CEDM and left-right four-quark operator, the expected dominance of nuclear effects comes from 
pion exchange due to both glr''^^ couplings, while only glr^ is present at LO for the ^-term. The isoscalar 
coupling glr^ approximately cancels in dh + dt, so while for the CEDM dh + dt ^ 3dd, for the ^-term 
dh + dt — O.S{dn + dp). Effects of the left-right four-quark operator can only be separated from CEDM 
if its TVPV three-nucleon potential is significant. For the quark EDM, where nuclear effects are much 
smaller, dh + dt — 0.8((i„ + dp) also holds but in addition one expects d^ — dt ~ 0.9((i„ — dp). The 
situation is most complicated for chiral-invariant sources, for which nuclear effects are significant for 
both A = 2,3, but they depend in the deuteron only on gl^^ while in the tri-nucleon g^ '' and Ci^2 
contribute as well. In this case dh + dt ~ 3dd — 2{dn + dp). By confronting these relations, measurements 
of light nuclear EDMs, particularly if they include the triton, could shed light on the mechanism of 
CPV [TU5] . 



4.2 Heavy nuclei, Shielding, and SchifF moments 

Nuclear physics is important in determining the EDMs of neutral atoms. And the primary fact from 
which all other considerations stem is expressed by the Schiff theorem |121j . which states that in the limit 
that in the limit of a point-like nucleus and non-relativistic electrons any nuclear EDM is completely 
screened by the atomic electrons, so that the net atomic EDM is zero. 

We give a brief illustration of this result, following Ref. xxx. Consider a system of structureless 
components (a nucleus and electrons), the k^^ of which has dipole moment dk, interacting via the 
Coulomb force V{r), so that 

k k k 

= Ho + i'^i^l/ek) ^dk-Pk,HQ^ ■ 

k 
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The perturbing Hamihonian (the last term above) shifts the unperturbed ground state |0) to 



\o) = \o) + J2 



|m) (m| Hd |0) 



En — Er, 



|o>+E 



|m) {m\iY,kO-/ek)dk ■ Pk |0) (-^o - -^n 
Eo — Em 



I l + z^(l/efe)d" 

V k 



k-Pk |0) 



(4.161) 



The induced dipole moment d' is then 

d'={0\J2^i^j\~0) 



z(0| 
-d. 



* ■ Pk 



|0) 



Y,dk 



(4.162) 



so that the net dipole moment of the entire system vanishes. The assumptions underlying this result 
are that the constituents are point-like, non-relativistic, and non-interacting except via the Coulomb 
force. In real systems, none of these assumptions hold fully. As we shall see immediately below, 
the finite nuclear size essentially leads to the replacement the nuclear dipole operator by the nuclear 
"Schiff operator," which contains two extra powers of the nucleon coordinate. Moments due to finite 
nuclear size are thus generically smaller by O (Rl^^i/Rl^^^) than the unscreened nuclear EDM. In 
diamagnetic atoms, the nuclear physics of which is discussed next, this suppression is mitigated by 
relativistic electrons and can be further mitigated by nuclear octupole deformation. In paramagnetic 
atoms, discussed in the next section, relativistic electrons can lead to a large enhancement of the atomic 
EDM. 

Further analysis leads to the result that the post-screening CP- violating nucleus-electron interaction 

is 



H = AnS- V6^r) + . . 



(4.163) 



where the omitted terms come from higher multipoles, e.g. the nuclear magnetic quadrupole (M2) and 
electric octupole (E3) multipoles. The operator S is the nuclear Schiff operator, defined as 



S = S'^ + s 



with 



'ch 



e 
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(4.164) 

(4.165) 



S 



N 



6 ^ 



(4.166) 



Here S'^^ is due to the charge distribution of the nucleus (usually the dominant piece), is due to 
the EDM of the nucleon, e is the charge of the proton, (r^)^^^ is the mean squared radius of the nuclear 
charge distribution, and dj is the EDM of nucleon j. The sum in Eq. (4.166) is over all nucleons, while 
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that in Eq. (4.165) is restricted to protons. Rotational symmetry lets us express the ground-state 
matrix elements of the three vector Schiff operators in terms of a single quantity: 



5 = (*o|5,|^o) , (4.167) 

where |\E'o) is the member of the ground-state multiplet with = J. 

The charge-distribution part of the Schiff moment, S'^^, can only be induced by an effective T— and 
P-violating inter-nucleon interaction. Most studies have been dedicated to the OPE part of the TVPV 



potential (4.158). The moment can have many sources, as we have seen, and can depend on other 
quantities besides the cjn^ . 



Equation ( 4.165[ ) is, as mentioned, only approximate. Corrections come from nuclear quadrupole 



deformation (which introduces a term proportional to the nuclear quadrupole moment), from relativity 
in electronic wave functions (which gives terms of order (Za)'^) |122l 1123] . and more subtle electron- 



nucleus interactions |124j . the complete forms of which are still not entirely settled. Equation (4.158) 
is also only approximate, representing the leading-order part of the chiral effective potential. Contact 
terms and higher-order pieces in effective field theory (which in heavy systems would be hard to control) 
or heavier-meson exchange in older frameworks will modify Vtvpv- At present, however, nuclear- 
structure theorists have not incorporated any of these corrections save (occasionally) those of order 
{ZaY into their calculations of Schiff moments. 

Beyond-the-standard-model and hadronic physics, as we have seen, determine the glr^ and the nu- 
cleon EDMs. The job of nuclear-structure theory, within the framework just defined, is to determine 
the dependence of the Schiff moment on these quantities. (Atomic physics in turn determines the de- 
pendence of the atomic EDM on the Schiff moment.) Here we examine only the dependence on the glr^ 
and the dependence on the nucleon EDMs can be computed as well, but is weaker. Only a few of 
the calculations cited below (e.g., Ref. |125j ) considers this weak dependence. We can parameterize the 
dependence on the as follows: 

S = (ao -g^:^ + a, g-g^:^ + a, g^^^) . (4.168) 

All nuclear structure information is thus encoded in the coefficients Oj, which have units efm'^. 

In what follows we discuss attempts to calculate the in several important nuclei. Most take 



advantage of the weakness of Vtvpv compared to nuclear energies and approximate S in Eq. (4.167 
essentially perfectly, by 

where |$o) is the "unperturbed" ground state — that obtained with Vtvpv turned off — and the 
are the corresponding excited nuclear states. 

4.2.1 i99jjg 

The atom associated with this nucleus has for years had the best limit on its EDM, and so ^^^Hg has 
received more attention by nuclear-structure theorists than any other nucleus (though still not nearly 
enough, as we argue below). Calculations range from the extremely schematic to the very sophisticated. 
The table below quotes the results of four, with brief (and inadequate) phrases signifying the techniques 
they employ. (A more extensive table, reporting several of the different estimates in, e.g., Ref. |125j 
as well as earlier versions of the ^^^Ra numbers presented in a later table can be found in Ref. 



The first nontrivial calculation was that of Ref. |127j : it approximated the unperturbed states in Eq. 



(4.169) by the eigenstates of a simple one-body potential and then treated Vtvpv approximately as 
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a zero-range interaction between the valence nucleon and the Hg core in first-order perturbation 
theory theory. Ref. |41| improved this treatment considerably by using the full finite-range Vtvpv 
and adding to the perturbative treatment the collective excitation of the core in the random phase 
approximation (RPA) by a simplified version of the residual strong nucleon-nucleon interaction. The 
resulting "core polarization" decreased the sensitivity of S to Vtvpv, as the table shows. Ref. [39] also 
used RPA to treat core polarization but in a diagrammatic version of self-consistent Skyrme mean-field 
theory (also known nowadays as energy-density-functional theory). The calculation, which contained 
of a self-consistent mean-field calculation in ^^^Hg before the treatment of core polarization, employed 
several state-of-the-art Skyrme energy-density functionals, giving rise to the range of numbers in the 
table. Finally, Ref. jl25] carried out the self-consistent mean-field theory (again with a number of 
Skyrme functionals) directly in the odd nucleus ^^^Hg, implicitly including the effects of RPA core 
polarization by the valence nucleon. It also allowed for axially-symmetric nuclear deformation and 
included O ((Za)^) corrections to the Schiff moment. Its result for the coefficient Oi is noticeably 
different from those of the other calculations, a fact that is hard to understand because the methods 
appear to include much of the same physics, albeit in quite different ways. 

Which of the calculations is most reliable and what is the uncertainty in our knowledge of the 
coefficients a,? Even if all the calculations included the same kinds of corrections to the naive Schiff 



operator in Eq. (4.165), these questions would be hard to answer. The calculations agree, more or 
less, on the size of oq and a2, but do not even agree on the sign of ai. Some possible reasons for the 
discrepancy between Refs. [39] and |125j . which, as mentioned, seem to include essentially the same 
many-body effects: 

(a) One of the calculations is in error. Ref. |125j carried out several internal consistency tests, but 
did not agree with Ph.D. thesis leading to Ref. [39] when repeating one of the calculations there. 
The authors suggest as a result that [39] may contain an error. On the other hand, the results of 
[39] agree fairly well with those of the similar RPA calculations in Ref. [H], suggesting that it is 
Ref. [125] , if any, that has problems. 

(b) Some of the mean-field solutions in Ref. |125j are metastable, though they are supposed to rep- 
resent stable ground states. But even those that are completely stable use the same Skyrme 
functional as [39], and — like the solutions in that reference — correspond to spherical shapes, 
disagree with [39] . 

(c) The treatment of core polarization in the two kinds of calculations are equivalent only up to 
terms linear in the strong interaction between the valence nucleon and the core, and only if that 
interaction is not density dependent. (It is in fact density dependent in Skyrme functionals.) But 
it is hard to imagine higher-order effects or the density dependence being very important. There 
are a few diagrams in Ref. [39] that have no counterpart in the odd-nucleus mean-field calculation 
of [125j . but their contributions are apparently small. 

(d) The state of the valence nucleon is represented only approximately in [3S]. Again, though, the 
approximation should be reasonably good. 

In short, it is difficult to see how the calculations could disagree so seriously. The authors need to revisit 
their work. 

It may, however, that all the calculations just reviewed are missing something important, and that 
the spread in results reflects their inadequacy. ^^^Hg is a soft nucleus in which a single mean field, the 
dominance of which underlies all of the results obtained thus far, is probably insufficient. It is, perhaps, 
unfortunate that ^^^Hg is such a difficult nucleus, but a better job is not beyond the means of nuclear 
structure theorists; techniques to mix mean fields with different properties exist. We contrast the state 
of affairs here here with that in another important and complex nuclear-structure problem: calculating 
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the matrix elements that govern neutrinoless double beta decay in complicated nuclei such as ^^Ge. 
Theorists believe that they know those matrix elements to within a factor of two or three, mainly by 
dint of the number of varied and careful calculations that have been carried out. The main theoretical 
problem with the Schiff moment in ^^^Hg is not the challenging nature of the calculation, but rather 
that only a few groups have tried. 

^^^Xe presents many of the same problems as ^^^Hg. Fortunately other nuclei, including the one we 
discuss next, are better behaved. 



4.2.2 



225 



Ra 



This nucleus became the focus of experimental interest after it was shown |128l 11291 1130] that the 
Schiff moments in nuclei with asymmetric shapes could be enhanced by two to three orders of magni- 
tude. ^^^Ra is octupole deformed, has favorable atomic physics, and has nuclear spin 1/2, making the 
nuclear orientation insensitive to stray quadrupole fields. It is thus presents a terrific opportunity to 
experimentalists. 

Shape asymmetry implies parity doubling (see e.g. Ref. |131j ). i.e. the existence in ^^^Ra of a very 
low-energy |l/2~) state (55keV |132j above the |$o) = 11/2'*') ground state, according to measurements). 



That low-lying excited state dominates the sum in Eq. (4.169) because of the small energy denominator 



it introduces. In the (good) approximation that the shape deformation is rigid, the ground state and 
its negative-parity partner in are projections onto good parity and angular momentum of the same 
parity-mixed and deformed "intrinsic state," which represents the wave function of the nucleus in its 



own body-fixed frame. Equation (4.169) reduces in these circumstances to |129] 



(4.170) 



c _ _2 ^ {Vtvpv) 
3^^^^55keV)' 

where the brackets indicate expectation values in the intrinsic state. 

The results of a couple of Schiff-moment calculations appear in Tab. 
|127] in ^^^Hg, obtained the intrinsic state by filling single-particle levels in a phenomenological octupole- 



10 Ref [123, much hke Ref. 



deformed potential and using a zero-range approximation to Vtvpv ^ but using Eq. (4.170) instead of 
summing over many unperturbed states as in Eq. (4.169). Ref. |10], like Ref. |125] in ^^^Hg, treated 



the (octupole-deformed) potential in completely self-consistent Skyrme mean-field theory with several 
Skyrme functionals, leading to a range of values for the Oj. It also included (and perhaps exaggerated) 
the damping effects of short-range nucleon-nucleon repulsion. Even so, the octupole deformation makes 
the resulting coefficients much larger than in Hg. 

As we have already mentioned, these calculations in Ra are almost certainly more reliable than 
those in Hg. The low l/2~ energy implies that the octupole deformation is strong and rigid, so that 
a single mean-field shape accurately represents the intrinsic density. There is thus little need to go 
far beyond mean-field theory here. Furthermore, experiments promise to increase the calculations' 
reliability. Theoretical work in progress jl33j shows that intrinsic Schiff moments are strongly correlated 
with E3 transitions, which have been measured in ^^^Ra [r34j at ISOLD E and m ay be measured in ^^^Ra 



itself. The resulting data will tightly constrain the factor {S^) in Eq. (4.170), leaving (Vrypy) as the 
only real unknown. Although a reliable calculation of that quantity is not trivial, it is far easier than 
calculating the transition matrix elements of both Vtvpv and 5*2 to all the excited states of ^^^Hg. 



4.2.3 Other Nuclei 

Theorists have calculated the Schiff moments of other nuclei as well, though not with as much care as 
they have in the nuclei already discussed. ^^^Xe has and will be the subject of experiments, and so has 
received some attention; like ^^^Hg, however, it is unfortunately soft. Researchers have also examined 
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actinides other than ^^^Ra, including some which have no static octupole deformation; the idea there is 
that dynamic deformation, i.e. octupole vibrations, may enhance the Schiff moments |135j . They have 
also considered the spherical nucleus ^^^Rn, which is to be examined experimentally as the first step in 



a project to work with heavier octupole-deformed Rn isotopes. Table 11 lists some of these results. We 



have omitted nuclei that show little prospect of being studied experimentally. 
4.2.4 Ranges and Best Values for Schiff Moment 



Table [12] lists best values and ranges for three important nuclei. We determined these in a somewhat 
subjective manner, assessing the strengths and weaknesses of each calculation. The entries should be 
considered tentative, and we cannot assign a quantitative meaning to our ranges; we simply consider 
it likely that the true values lie in them. For the case of ai in ^^^Hg, as already discussed, our range 
includes zero. 



5 CP and T at the Atomic and Molecular Scale 

As with the physics at the hadronic and nuclear scales, it is convenient to express the atomic and 
molecular EDMs in terms of the operators that characterize physics at shorter distance scales. To 
that end, we first write down an expression for a general atomic or molecular EDM (Ia in terms of the 
electron and nucleon EDMs, the nuclear Schiff moment, and the Wilson coefficients for the dimension 
six T- and P-odd electron-quark interactions. We subsequently express the electron EDM in terms of 
either the Wilson coefficient or the quantity 5e- Doing so allows us to express dA in such a way 
as to place all of the fundamental dimension six operators on the same footing, and in the case of the 
electron EDM, take into account the additional Yukawa suppression that accompanies this operator. 
Thus, we have 



dA = p\de+ pzdN + KsS + 

N=p,n 



A 



+ 



4%f + k?g? 



k?g? + k^^^g? 



ImC(-) 



(5.171) 
(5.172) 



where we may alternately express the electron EDM contribution as 



(5.173) 



Note that the d^ and Schiff moment S may then be expressed in terms of 6 and the dimension six quark 



and gluon operator coefficients using Eqs. (3.42 3.46) and (4.168), allowing one to explicitly identify 
a common factor of (w/A)^ for all dimension six operators and, thereby, to place them on a similar 
footing. 

To illustrate the relative sensitivities of various atomic and molecular EDMs to the quantities appear- 
ing in Eqs. (5.171 5.173) we consider one paramagnetic atom (^°^T1), one diamagnetic atom (^^^Hg), 



and one molecule (YbF) for which the most stringent experimental limits have, thus far, been obtained. 
A summary for other cases is given in Table. As a prelude, we first summarize a few features of the 
atomic and molecular computations, referring the reader to the extensive reviews in Refs. [TT| [T3] and 
the recent study in Ref. |136j for details. 

The sensitivity of an atom of molecule to the electron EDM is governed by corrections to the Schiff 
screening as describe above . In contrast to the corrections due to finite nuclear size or higher T- and 
P-odd nuclear moments, the relevant corrections are relativistic and entail both a first and a second 
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order energy shift proportional to de- ^E^f) for J = 1; 2. For our purposes, the expHcit expressions are 

not particularly enlightening, and we again refer to Refs. [131 1124] for details. The correction Ai?|^j 
can be expressed alternately in terms of an electronic operator proportional to 75 or 70 — 1, indicating 
the dependence on lower components of the electronic wavefunction that vanish in the non-relativistic 
limit. The second order term ^E^^-^ entails polarization of the atomic cloud by the presence of de (again 
vanishing in the non-relativistic limit) that is then probed by the external field. For heavy paramagnetic 
atoms, the polarization correction dominates, growing as Z^. 



The four-fermion, semileptonic interactions in Eq. (3.40) lead to an effective atomic Hamiltonian 



that takes on the following form in the limit of an infinitely heavy nucleus: 



fVatom 
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(5.174) 

(5.175) 
(5.176) 

(5.177) 



and where the Dirac matrices, 5{r) and V act on the electronic wavef unctions. 

Note that in arriving at the expression for Hs we have performed the sum over all nucleons, using 
the fact that in the non-relativistic limit the operator NN just counts the number of nucleons. For 
Ht,p, in contrast, the nuclear matrix elements of the spin operators is more complicated. The results 
quoted below for heavy nuclei, which have N > Z, were obtained assuming a single unpaired neutron 
contributes and using a single particle shell model result for the nuclear matrix element of Thus, 
the values for the A;^p quoted below correspond only to the neutron contribution (or the difference 



k. 



(0) 



k 



(1) 



j^p — rvrppj. For all paramagnetic atoms, all three interactions Hs^p^t contribute. For diamagnetic 

atoms wherein all electrons are paired, Hs cannot induce an EDM except in tandem with the hj^erfine 
interaction. 



5.1 Paramagnetic atoms: Thallium 

According to the computations of Refs. [HI [13] (see also Ref. |136j for a recent summary), the EDM 
of 205x1 has by far the strongest dependence on the electron EDM and the ei'y5eNN interaction of all 
paramagnetic atoms studied experimentally to date. Compared to the latter, the dependence on tensor 
and nucleon pseudoscalar operators are suppressed by three and five orders of magnitude, respectively. 
Numerically, one has |137] 

p^ = -573±20 , /37 = 0.65 ± 0.02 e fm , = (1.9 ± 0.07) x 10"^ e fm , (5.178) 

while 

kP = -(7 ± 0.3) X 10-^ e fm and = 0-24°^ , (5.179) 

and JT3\ 

kP = -kP = 1.5 X 10"^° e fm and k^ = -k^ = 0.5 x 10"^ e fm . (5.180) 

The numerical dominance of kf^ implies that (iA(^°^Tl) has the greatest sensitivity to lmCeq \ a some- 
what reduced sensitivity to ImCi^'* and relatively little sensitivity to ImC^^g]^. 



45 



It is interesting to compare the relative sensitivity of (iA(^°^Tl) to 5e and ImC, 



(-). 



eq 



, (0) 

^^-37 . (5.181) 

To the extent that these two quantities have the same order of magnitude, the four-fermion semileptonic 
operator would yield a far larger contribution to the thallium EDM than would the electron EDM. 
The corresponding sensitivities for i33Cs[l38l |l39l [HQ], ^^Rb [I39], and ^lOprlHTl Ell [Hal [I36] are 



also listed in Table 13 We have largely followed Ref. |136j in averaging the results for Cesium and in 



assigning error bars to the results for Francium, for which only the analytic expressions in Ref. |142j 



have been used to obtain k!f\ 
5.2 Mercury 

The diamagnetic nature of ^^^Hg makes it far more sensitive to T- and P-odd interactions at the purely 
hadronic level than those involving electrons as compared to the paramagnetic thallium atom. From 
Ref. Refs. [TT|, [T3] we first obtain the sensitivity to the nuclear Schiff moment: 

^5 = 2.8 X 10"^ fm-2 (5.182) 

where the scale for the Schiff moment is e-fm^. While it is possible to include the nucleon EDM contri- 
butions in the nuclear Schiff moment, we find it helpful to separate these contributions out explicitly 



as in Eq. (5.171). From Ref. pT| [13] we obtain 

pP = -5.6 X 10"^ p\ = -5.32 X 10"^ , (5.183) 

while the sensitivity to the electron EDM is given by 

= 0.01 , /37 = -1.13 X 10"^ e fm , €(% = -3.3 x 10"" e fm , (5.184) 

essentially five orders of magnitude less sensitive that 205'ji_ Tj-^g sensitivity to the four-fermion semilep- 
tonic interactions NNei'j^e are similarly suppressed with respect to thallium: 

kPCRg] ^ 1.16 X 10-^ X fcf[205Tl] , (5.185) 

while the sensitivity to the tensor and nucleon pseudoscalar interactions are somewhat enhanced: 

fcg^[i99Hg] ^ 4xfcg^[205Tl] , J = 0,1 (5.186) 
4')[^99Hg] ^ 4x4'^p05Tl] , J = 0,1 . (5.187) 

The relatively weak dependence of (i^(^^^Hg) on the scalar interactions reflects the suppression due to 
the presence of the atomic hyperfine interaction that must be present to yield a no n- vanishing result. 
Consequently, (iA(^^^Hg) provides a relatively more effective probe of ImC^f^^ than does (i^(^°^Tl). In a 
scenario where the only T- and P-odd effects arise via semileptonic interactions, a comparison of results 
from mercury and thallium could allow one to disentangle between various sources. 

Before proceeding with the molecular EDMs, we observe that if a given BSM scenario generated 
only the dimension six quark EDM operators and not the four-quark, CEDM, or three gluon operators, 
then the corresponding effect on ^^(^^^Hg) would be dominated by the induced neutron and proton 
EDMs. In this case, the present limit on (iyi(^^^Hg) could be interpreted as a limit on dn at the 10~^^ 
e fm level, roughly one order of magnitude weaker than the present direct neutron EDM limit. On the 
other hand, using the latter, one could then infer a bound of roughly 10~^° — 10~^^e fm on dp. 
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Figure 4: Dependence of ri{Eex_t) for YbF |145j . Vertical axis gives E^s as a function of the 
applied field E^xt- Scaling E^g by twice its asymptotic value gives rj{Ecxt) ■ Figure reprinted 
with permission from Phys. Rev. Lett. 89 023003 (2002). 



5.3 Ytterbium Fluoride 

In order to assess the sensitivity of the polar molecules to the underlying CPV operators, we first 
convert to the conventions used in the theoretical literature. Following Ref. |144j we write the molecular 
Hamiltonian as 

HJ^oi^y = {Wde + WsCs) S ■ h (5.188) 

— * — 

where S is the electronic spin, n is a unit vector along the axis of the YbF molecule, and Cs = 
{Z + N)Cf^ + {Z — N)C^\ We do not at present include the dependence on the pseudoscalar or tensor 
interactions as, to our knowledge, the corresponding evaluations of molecular sensitivities have not been 
performed. 



Contrary to what one might naively expect, the interaction in Eq. (5.188) contains no dependence 
on the external electric field, -Eext. The experimental observable - a frequency shift - depends on the 
ground state (g.s.) expectation value in the presence of the external field: 

(g.s.l H^^'py |g.s.)^^^^ = {W,d, + WsCs) r/(Eext) , (5.189) 

where 

v{E,xt) = {g-s.\S-n\g.s.)j,^^^ (5.190) 
is an effective polarization that increases monotonically with E^xt = \Eext\ and has a maximum value 



of 1/2. From the first term in Eq. (5.188), then, one may interpret ri{Ecxt)Wd as the effective internal 
electric field E^s acting on an unpaired electron that is induced by a non- vanishing -Eext- For YbF, 
ri{Eext) has been reported in Ref. [145J and is shown in Fig. 111. The vertical axis gives E^s as a function 



12 



of Eext- The value of ri{Eext) may be obtained by scaling E^s by twice its maximum, asymptotic value 
For purposes of this review, it is useful to express H^ypy in terms of 6e and the ImCef^: 

HTVPV={iy [eCASe + gfh^hmCi;^ + g'ih'ihmCip] S-h , (5.191) 

where the quantities CI and k^g'^^ are determined by molecular structure. The results for and Ws 
are typically quoted in units of Hz/(e cm) and cm, respectively . For the sensitivity to the electron 
EDM, the latest results in the compilation of Ref. |136] yields 

l^rf = -(1.1 ±0.1) X 10^2 Hz/(e fm) , /3^^ = (1.2 ± 0.1) x 10^ Hz , eCl = 3.6 ± 0.3 kHz (5.192) 



^We thank T. Chupp for a helpful discussion of this point. 
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where we identify p\ = Wd- For the scalar interactions one has 

kP = 54^^ = -(92 ± 9) kHz , (5.193) 

where we have adopted the ten percent theoretical error suggested in Ref. |136] . To compare with 
the thallium atom, one has ^V^Cl ~ —26, indicating a somewhat stronger relative sensitivity to 
the electron EDM. However, the sensitivities are sufficiently similar that a combination of the present 
experimental limits in the two systems does not allow for a significant individual limits on 6e and 

Looking to the future, an effort to probe the EDM of ThO is underway. A value for Wd has been 
computed in Ref. [146] . The corresponding value for Ws has been inferred from t 
computed analytically in Ref. |142] . A conversion to e^l and k^p appears in Tables 
cases, we have arbitrarily assigned a ten percent theoretical uncertainty. We note that the ratio k^g' /e(% 
for ThO lies approximately midway between that of Th and YbF. 
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6 Beyond the Standard Model: Examples 

The space of BSM scenarios that contain additional sources of CPV is vast, and it is not feasible 
to provide an exhaustive review here. Instead, we will focus on several representative examples to 
illustrate the interplay of scales: supersymmetry, models with extended gauge symmetry, and scenarios 
with extra spacetime dimensions. Before doing so, we first make a few general remarks. Perhaps most 
importantly, any new source of CPV will generally induce a contribution to the QCD vacuum angle, 
which we denote as ^bsm- In the minimal supersymmetric Standard Model (MSSM) for example, such 
contributions arise at one- loop order via corrections to the quark propagators. Given the already severe 
bounds on 9, such contributions to ^bsm by themselves imply stringent limits on the CPV phases in the 
absence of a mechanism to alleviate them. Possibilities include invoking a new symmetry, such as the 
PQ symmetry or a flavor symmetry that yields a vanishing one-loop result. 

Second, non-observation of atomic, molecular, and neutron EDMs generically imply that any new 
CPV phases 0cpv must be quite small if the BSM mass scale A is sub-TeV. Conversely, allowing 
sin^cpv ~ 1 implies that A > few TeV. To illustrate, dimensional analysis gives for the elementary 
fermion EDM 

^/ ~ ^ {ji) 2 sin^cpv (6.194) 
where ak is either the fine structure constant or strong coupling (evaluated at the scale A). For ak = aem 



Eq.( |6.194D gives 

-^j xlO-^=^efm . (6.195) 

The present limit on the EDM of the electron, \de\ < 10.5 x 10^^^ e fm |147] obtained from an experiment 
on the Yb-F molecule, then implies that 

|sin0cpv| < (A/2 TeV)^ . (6.196) 

Thus, for I sin0cpv| ~ 1 one requires A > 2 TeV. In order to allow for sub- TeV scale masses and 0{1) 
CPV phases while respecting present constraints, one must either invoke cancellations between different 
contributions Pl] or a mechanism that suppresses the one-loop EDMs. In the case of the MSSM, 
for example, taking the sfermions to have mutli-TeV masses can result in the leading contributions 
arising at two-loop order and involving the electroweak gaugino-Higgs/Higgsino sector with sub- TeV 
masses [EHl 1149] . Given the suppression of an additional loop factor, the resulting dependence on the 
CPV phases is weakened and the present constraints are generally less severe. Alternate strategies can 
be employed in other scenarios. 
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6.1 Supersymmetry 

Supersymmetry (SUSY) remains one of the most strongly motivated BSM scenarios, providing an 
elegant solution to the hierarchy problem, candidates for cold dark matter (the lightest neutralino or 
gravitino), and copious sources of CPV that may drive the generation of the baryon asymmetry during 
the EWSB era. At the same time, SUSY CPV generically leads to one-loop EDMs that exceed present 



experimental bounds, assuming that superpartner masses lie below one TeV [see Eq. (6.196)], leading to 
the so-called "SUSY CP problem". The SUSY mechanism for solving the hierarchy problem leads one 
to expect sub- TeV scale superpartner masses, implying | sin^cpvl ~ 0.01 — 0.1. On the other hand, one 
might naturally expect sin^cpv ~ C^(l)- Moreover, successful supersymmetric electroweak baryogenesis 
typically requires C(l) phases unless the relevant portion of the superpartner spectrum is finely-tuned 
to contain near degeneracies. 

Several solutions to the SUSY CP problem have been proposed: 

(i) Heavy sfermions. It is possible, for example, that the fermion superpartners (sfermions) are 
considerably heavier than one TeV, leading to a suppression of one-loop EDMs and allowing for 
C(l) phases p30| I148[ 1149] . Null results for superpartner searches at the LHC may be pointing to 
this "split SUSY" scenario, as the generic mass bounds on gluinos and first and second generation 



squarks are now at the TeV scak^ The electroweak gauge boson and Higgs boson superpartners 
may still be relatively light, thereby allowing for a viable baryogenesis mechanism (see Ref. [1] 
and references therein). 

(ii) CP- conserving SUSY breaking. It is equally possible that the mechanism of SUSY-breaking that is 
responsible for both the splitting of SM masses fro those of their superpartners (the "soft terms" ) 
and the CPV phases suppresses the latter. This possibility has been emphasized in the work 
of[151] that considered an M-theory scenario on the G2 manifold in which the only source of CPV 
at the SUSY-breaking scale is the CKM phase. The resulting effects on low-scale parameters then 
enters through the RG evolution. 

(iii) Cancellations. It was proposed some time agofS^ that contributions to EDMs from different CPV 
phases or those from different dimension-six CPV operators may cancel leading to a suppression 
that again allows for 0{1) phases and light superpartners. 

An extensive discussion of (i) and (iii) are given in the reviews of Refs. 126] and the more recent 
analysis of EDMs in SUSY given in Ref. |88]. Given the comprehensive nature of these articles, we do 
not attempt to provide an exhaustive review of EDMs in SUSY here. Instead, we summarize several 
generic features as well as developments that have appeared since publication of these studies. For 
this purpose, we focus on the minimal supersymmetric Standard Model (MSSM) for which one has the 
superpotential from which one derives the supersymmetric Lagrangian, 

W^MSSM = uyuQHu - dydQHd - eyeLHd + fJ-Hu ■ Hd- (6.197) 

Here, the hatted quantities / and F are the SU(2)2,-singlet and doublet chiral superfields for fermion 
F while are the two Higgs doublet superfields. The y/ are 3x3 Yukawa matrices. For purposes of 
this discussion, we omit possible R-parity violating terms in the superpotential, which is tantamount to 
promoting the accidental global B — L conservation of the SM to a symmetry of the MSSM. Note that 
superpotential introduces only one new parameter beyond that of the SM, namely, the coefficient of 



the last term in Eq. )(6.197). In addition, EWSB allows the two neutral Higgs scalars to have vacuum 



expectation values, whose ratio defines the angle (3: tan/3 = Vu/vd. 

^^li is still possible, however, that these strongly interacting superpartners are lighter than one TeV but have a 
compressed spectrum leading to presently undetectable experimental signatures. 
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The soft SUSY-breaking Lagrangian responsible for splitting the SM and superpartner masses is 

Aoft = -]^{M^r9 + M2WW + M^BB) + c.c. 

-{usluQHu - dsLdQHd - eaeLHd) + c.c. 

— Q^nigQ — Umf^L — Mm?M^ — Jm|(J^ — eriige^ — 7Tijj^H*Hu — rn^jj^H^Hd 
-{hH^Hd + c.c.) (6.198) 

Here, the first line gives the gaugino mass Mj, i = 1,2,3 for the U(l)y, SU(2)l and SU(3)c gauginos, 
respectively. The second line gives the trilinear "A-term" that couples Higgs scalars with left- and 
right- squarks and sleptons. The third line gives the scalar mass tn'-^ ^, rnj , and m\ ^ for squarks, 
sleptons and Higgs scalars, respectively. As with the Yukawa matrices, the boldfaced quantities indicate 
matrices in flavor space. Finally, the last line is the bilinear 6-term, which couples up- and down-type 
Higgs scalar doublets. It is also important to emphasize that in nearly all SUSY analyses, one takes 
the a/ to be proportional to the corresponding Yukawa couplings, thereby naturally suppressing flavor 
changing neutral currents: 

= Afjf (6.199) 

where Af becomes the effective trilinear soft parameter for each fermion species. 

The various interactions in £soft introduce copious sources of both flavor violation and CPV. Here, 
we focus on the latter. After performing an appropriate set of field redefinitions, £soft - together 
with the /x-term in the superpotential - includes 40 CP-violating phases beyond those of the SM (for 
a useful discussion, see, e.g., Ref. |152] ). Unlike the CPV phase in the CKM matrix, the effects of 
these new phases are not suppressed by the Jarlskog invariant [153J and light quark Yukawa couplings. 
Consequently, the CPV effects need not be suppressed as in the SM, leading to the SUSY CP problem. 
The new phases can be classified in terms of those that solely enter the gauge-Higgs sector: 

0, = Arg (/iMi6*) 0,, = Arg (M.M*) (6.200) 

where i, j run over the three gauge groups of the MSSM (leading to a total of three independent phases 
in this sector which we take to be the (pi ), and those involving the sfermions: 

0/ ^ Arg {AfM*) <Pff, ^ Arg {AfA},) . (6.201) 

Given the large number of phases, one often invokes a phase universality assumption: 

01 = 02 = 03 = 0M (6.202) 
0^ = 0^, = 0^ . (6.203) 

As we discuss below, this assumption is unlikely to allow for consistency between supersymmetric 
electroweak baryogenesis and EDM constraints. 

The CPV interactions in the MSSM give rise to three of the dimension six operators of interest in 
this article: the EDM and CEDM operators arise at the one-loop level, while the three-gluon operator 
first enters at two-loop order. The four- fermion operators are technically dimension eight, but can be 
enhanced for large tan /3. Illustrative contributions to the one loop EDM and CEDM arise from the 
diagrams in Fig. [Sj 

Although not shown explicitly, the insertion of the Higgs fields as needed for electroweak gauge 
invariance enters in one of two ways: (a) the mixing of left- and right-handed sfermions and (b) mixing 
of electroweak gauginos and Higgsinos. The former is characterized by the sfermion mass-squared 
matrix: 

- ( m|: ) (''•^«) 
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Figure 5: Illustrative one-loop (left) and two-loop (right) contributions to the fermion EDM 
and quark CEDM in the MSSM . 



with 

ML = niQ + mJ-FAf (6.205) 

M^K = m? + + Af (6.206) 

with 

Af = (^4 - Qf sin^ 9yv^ cos 2/3M| (6.207) 



and 



Af = Qf sin^ 9w cos 2/3M| (6.208) 
, ^ \ V [af sin d — uvf cos B] , u — type sfermion 

M^K = M^L = i . • « J ; f • • (6.209) 

I f [af cos p — /iyf sm pj , d — type siermion 

Here niq is the mass matrix for the corresponding fermion; // and Qf are the third component of isospm 

and fermion charge, respectively; and v = ^/f^ + f^. The explicit factor of v in the M^j^ = M^j^l 
corresponds to the Higgs insertion above the scale of EWSB and leads to mixing between superpartners 
of the left- and right-handed fermions after EWSB. As a result the incoming and outgoing fermions in 
Fig. [5] can have opposite handedness as needed for the EDM and CEDM operators. Note that with 



the assumption of Eq. (6.199) the left-right mixing is proportional to the fermion Yukawa coupling, 



implying that the contributions to the EDM and CEDM are as well. 
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A similar Higgs insertion is implicit in the mixing of the Higgsinos and electroweak gauginos. To 



illustrate we give the chargino mass matrix for the charged fields ijj = {W^ 




X 



XT 





X 



M2 

V2cbMw 



(6.210) 



whose mass eigenstates are the charginos xf^ i = 1,2. Note that since = gv/2, the off-diagonal 
terms responsible for gaugino-Higgsino mixing contain an implicit Higgs insertion. Thus, an incoming 
left-handed fermion in Figjsjthat interacts with the charged SU(2)i gaugino (the "wino") component of 
the xf can lead to an outgoing right-handed fermion that interacts with the Higgsino component due 
to this mixing. Since the latter interaction is given by the Yukawa interaction in the superpotential, the 
corresponding effect on the EDM and CEDM is again proportional to the fermion Yukawa coupling. 

A comprehensive set of expressions for the one-loop contributions to the fermion EDMs and quark 
CEDMs are given in Ref . [88] , so we do not reproduce them here. I nstead , we give an illustrative set of 
expressions for 6f and 6q under the universality assumption of Eq. (6.202): 



Qe 

327r2 

327r2 
1 

~32^ 



— sm (pA — 1 sm (p., 

12 V 24 24 ' 



24 ^ 24 ) si^*^/^^^^^ 
(sin0^[tan/3]±i - sin^^) + 0{glgl 

^ (sin0^[tan/3]±i - sin^^) + 0{glgl 

io 



(6.211) 
(6.212) 
(6.213) 



where Q/ is the fermion charge, where for simplicity we take SUSY mass parameters to be identical 
(A = \Mj\ = |/i| = I A/I), and where the upper (lower) sign corresponds to negatively (positively) 
charged quarks. 



The expressions in Eqs. (6.211) contain a linear combination of the two universal phases, allowing 
for the possibility of some cancellation between various contributions. However, as noted in Ref. [E^ 



it is unlikely that such a cancellation could allow one to evade all EDM limits since the coefficients of 
sin0^ and sin0^ differ for the various species of fermions as well as between the EDM and CEDMs. 
Nonetheless, it is in principle possible to obtain a consistent fit to present EDM limits with 0{1) sin 0cpv 
if one relaxes the universality assumption, a feature we discuss below. 

Going beyond one-loop order, one encounters the first contributions to the three-gluon operator 
as indicated in Fig. [6] as well as the two-loop Barr-Zee graph contributions to the EDM and CEDM 
operators as indicated by the diagrams in Fig. [sj^right). In the limit of heavy sfermions, the Barr-Zee 
graphs containing electroweak gauginos in the upper loop will give the dominant contribution. Explicit 
expressions for the contributions generated by the graphs of Fig. [6] are given in Ref. [88] , along with 
those for a subset of the Barr-Zee graphs that contain exchanges of only the lightest CP-even Higgs 
scalar. In Ref. [r49i| . the remaining set of graphs containing electroweak gauginos and exchanges of the 



charged and CP-odd Higgs as well as other gauge bosons were computed and found in some cases to 
give the dominant contributions to the fermion EDMs in the heavy-fermion regime. We note that in 
this regime, even the two-loop CEDM operators are suppressed, since the upper loops in Fig. [sj^right) 
contain only squarks. 

Based on that work, the authors subsequently performed a global analysis of EDM constraints on 
CPV phases in the MSSM [TM] . up-dating the SuperCPH2.0 code described in Ref. [SB] to include the 
full set of two-loop graphs. Illustrative results (obtained before publication of the YbF molecular EDM 
result) are indicated in Fig. [T] and Table 14 based on the use of QCD SR to compute the hadronic 
matrix elements. In this context, the impact of the three-gluon operators is typically suppressed, as is 
the effect of the four-fermion operators in the low-to-moderate tan (3 regime. The results in Table 14 
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Figure 6: Illustrative MSSM contributions to the CPV tliree-gluon operator 



were obtained assuming three independent phases contribute: 0i, 03, and a common triscalar phase for 
the first generation squarks: 0^ = 0^. The impacts of 0e and 0i are sufficiently weak that one may 
omit them from the global analysis, though the sensitivities of on-going and future EDM searches could 
allow one to probe these phases as well. Fig. [7] shows the relative correlations between pairs of phases, 
obtained from a fit in each case including only those two phases. 

It is particularly notable that the (iA(^^^Hg) limit places severe constraints on 03 while generating 
a strong correlation between this phase and (f)u,d, both of which enter the CEDM operators at one-loop 
order. In contrast, the neutron and Thallium EDM limits have a relatively stronger impact on 02, 
though at present the latter constraint is not strongly correlated with any of the other phases. Future 
measurements with ~ 100 times better sensitivity, however, would give rise to such correlations. 

6.2 Extended Gauge Symmetry 

As with SUSY, the embedding of the SM gauge symmetry in a larger gauge group can allow for additional 
CPV phases in both flavor diagonal and flavor non-diagonal processes at low energies. For purposes 
of illustration, we consider the well-studied left-right symmetry model (LRSM) with the gauge group 
SU(2)LxSU(2)ijxU(l)B_L. Symmetry breaking proceeds in two steps, with the first step breaking the 
left-right symmetry and generating a mass for the right-handed gauge bosons, followed by a second step 
that breaks the SM gauge symmetry. Implementing this scenario requires augmenting the SM Higgs 
sector with additional scalar fields: two complex triplets A^ r that transform separately under the 
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Figure 7: Constraints on MSSM CPV phases implied by null results for the neutron, Thal- 
lium, and Mercury EDMs flM] . Red, blue, and black contours correspond to first generation 
sfermion masses (M^^ /j/^)! = 200, 500, and 1000 GeV, respectively. With kind permission 
from Springer Science+Business Media: J. High Energy Physics, "A comprehensive analysis 
of electric dipole moment constraints on CP- violating phases in the MSSM", 08, 2010, p. 
062, Y. Li, S. Profumo, M.J. Ramsey-Musolf , Fig. 10 (partial). 



SU(2)x, /J symmetries and an eight-component bidoublet 0. New CPV phases arise from two sources. 
The extended gauge symmetry allows for a complex phase a associated with the VEV of </>, corresponding 
to spontaneous CPV (SCPV): 

= ( ) ■ 

In addition, new phases can arise in the fermion mixing sector owing to differences between separate 
rotations of the left- and right-handed fermion fields, Im(V^^V^^*). Significantly, one requires only two 
generations of fermions rather than three as in the case of the SM in order to obtain a CPV phase 
associated with mixing that cannot be removed through field redefinitions. 

The manifestation of these new phases then enters through the mixing of the Wl^h gauge bosons 
into the mass eigenstates Wi^2'- 

Wt = cos^W^ + sm^e-'''W^ (6.215) 
W+ = -sinCe^°Vr+ + cosCl^+ (6.216) 



where the mixing angle is given by 



tane = -«:«:74^-2- f ^ ) (6.217) 



with f /J being the vev of the neutral component of the Aji and being the mass of the boson. 
Interactions of the latter with quarks and leptons can yield a variety of the dimension six effective 
operators introduced earlier. The resulting contributions to the dg were first computed in Ref. |155] for 
the two-flavor case. In contrast to the SM, non-vanishing contributions first arise at one-loop order, 
with the result 

Sq = (^1^) sin ^ Ag [rq cos 9l cos 9r sin a + rq sin 9l sin 9r Bg] (6.218) 
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where cos Ol^r is the "Cabibbo angle" for the left- and right-handed sectors, h are the corresponding 
gauge couplings for the two sectors, 



ru = YjYu rQ = YjY^ A^ = A B, = sin{a + 5r- 5l) (6.219) 

rd = YjYd rQ = YjYd = 5 B, = sin{a - 5r + h) (6.220) 

and where we replace 

in the definition of 5q. Since Mi is the mass of the lightest W^-boson, one encounters no explicit 
suppression due to the heavy scale. Note also that the EDM for a given quark flavor is proportional to 
the Yukawa coupling for the quarks having opposite sign third component of Il,r- 

The corresponding CEDM operator has been computed in Ref. |156] One also encounters the four- 



quark operator of Eq. ( 2.29 ) due to the exchange of the Wi^2 between quarks |157t I156i 1271 128] . Following 



the notation of Ref. ([28j) one obtains 

^' ^ GpK^'^ sin^Im (e-'^V^^/^*) (6.222) 



A2 



where we have included the contribution from only the exchange of the Wi and have extended the 
operator Q^ud to include all down-type quarks q = d,s, b. The constant K^~^ ^ 3.5 is a QCD renormal- 
ization group factor associated with running from the weak scale to the hadronic scale. Again specifying 
to the two generation case, we observe that the effect of the phase in the quark mixing matrix will not 
enter the operators containing only u and d quarks, leaving only a dependence on the SCPV phase. In 
this case, the contribution to gi^^ will depend solely on this phase and not on 6l — Sr. Including the 



second generation quarks would then require extending the arguments leading to Eq. (3.151 ) to account 
for the nucleon matrix element of the ss and a coupling of the nucleon to the rj meson. We leave this 
extension, as well as a consideration of the CEDM and three-gluon operators, to future work. 

Contributions to the neutron EDM in the LRSM have been carried out using a variety of approaches. 
Ref. |155] relied on the quark model result to determine the dependence of dn on the dq. The authors 
of Ref. |157] also performed a quark model evaluation of the contribution from the four quark operator 



(2.29). Pseudoscalar loops were included in Refs. [156|(28], where one of the pseudoscalar meson-baryon 
vertices are induced by the underlying CP-violating quark and gluon operators, while Ref. also 
computed contributions to the nucleon wavefunction due to the CEDM. We note that the pseudoscalar 
loop results in Refs. |156l [28] were not performed using a consistent chiral power counting and are likely 
to overestimate the corresponding contribution that is proportional to 

To illustrate the manifestation of LRSM CPV in EDMs, we consider (a) contributions from the dq 



to dn using the quark model relation 3.138 (b) chiral loop contributions to dn induced by Q^uq', (c) 
contributions to (i^(^^^Hg) generated by Q^uq via the nuclear Schiff moment. Starting with the dq, we 
neglect the heavy quark contributions for simplicity and take Ql = gR = e/ sin6'w, leading to 

dn ~ (1.13 X IQ-^e fm) M - j sin^ (5r„ + 4^^) cos cos sin a . (6.223) 

Noting that bY^ + AUd ~ 2.5 x 10"^ and that | sin^| < 10"^ from tests of first row CKM unitarity[ISS], 
we see that 

\dnf^ < (3 X IQ-^^e fm) M - j^j cos cos sin a . (6.224) 
The present dn constraint is, thus, not sufficiently stringent to probe this contribution. 
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A potentially larger contribution may arise from chiral loops involving the isovector TVPV ttNN 
interaction. Making the same simplifying assumptions used above and using Eqs. ( |3l5l1[6:222| one has 

~ _10-4 h - — Lj sin^ cos^icos^ijsina . (6.225) 

The corresponding one-loop contribution to dn is given by 

where fin = —1.91 is the neutron anomalous magnetic moment and F{x) is a loop function. An early 
calculation reported in Ref. |156j gave F{x) = 3/2 — x + ■ ■ ■ , where the "+ ■ ■ ■ denote contributions 
non-analytic in x. This computation, however, did not utilize the consistent power counting obtained 
with HBxPT and, thus, should be considered unreliable. A consistent HBxPT computation gives 
F{x) = —xlnx ~ —0.1, implying an order of magnitude smaller neutron EDM contribution than one 
would infer from the computation of Ref. |156] . Taking | sin^| < 10^^ we then obtain 

1- cos^LCos^ijsina , (6.227) 

indicating roughly four orders of magnitude greater sensitivity to sin a than implied by the quark 
EDM contribution. One may trace this difference to the combination of the quark Yukawa couplings 
5Yu + 4:Ud ~ 2.5 X 10~^ that enters the quark EDM contribution and that does not appear in the 
chiral loop contribution induced by Q^^uq- For cos 9l ~ cos9r ^ 1 and Mi << M2 we then obtain 
I sina| < 10~^ from this contribution. 

Turning to ^^(^^^Hg), we use the value of ks given in Eq. (5.185) and a representative value for ai 
of 0.03 (midpoint of the corresponding range in Table 12) to obtain 

\dAC^^B.g)\ < (1.1 X 10""e fm) M - j^j cos cos ^j? sin a , (6.228) 

giving an even stronger sensitivity to the SCPV phase than dn, though subject to considerable nuclear 
theory uncertainties associated with the computation of ai as discussed above. However, given that the 
current bound on dAC^^Hg) is three orders of magnitude smaller than the limit on dn, the former is 
likely to provide the most stringent constraint on the LRSM contribution even allowing for a possibly 
smaller magnitude for oi than assumed here. For the benchmark value of ai used in this example, we 
would obtain | sina| < 10~^ for cosOl ~ cosOji ~ 1 and Mi << M2. 

In the foregoing discussion, we have used the phenomenological constraint on the mixing angle ^ 
obtained from tests of first row CKM unitarity |158j . An alternate approach has been followed by the 
authors of Ref. [28], who observed that one may determine the elements of the right-handed CKM 
matrix Vj^ in terms of Vj^^, the ratio k'/k, sin a and the quark masses by exploiting properties of the 
LRSM Yukawa matrices, the hierarchy of quark masses, and the Wolfenstein parameterization of V^^. 
One then finds that 

sineim(e-^"yiV;f) (6-229) 
can be expressed in terms of M^/M| and rsina, where r = {mt/mb){K,' / k) characterizes the ratio of 
the two bi-doublet vevs. Illustrative constraints on M2 ~ and r sin a are indicated by the yellow 
points in Fig.js] We note that the yellow points were obtained using the value of F{x) given in Ref. |156j 
that is an order of magnitude larger than the HBxPT result. Taking into account the latter reduction 
and utilizing the bounds on dAC^^Hg) we conclude that the region allowed by the ^^^Hg limit is likely 
to be considerably narrower than indicated by the yellow points in Fig. [8j 

A future improvement in the sensitivity of neutron EDM searches by two orders of magnitude could 
make dn a comparably powerful probe of LRSM CPV as ^^(^^^Hg). 
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Figure 8: Constraints on right-handed VT-boson mass and CPV parameter r sin a obtained 
from present constraints on the neutron EDM (yellow dots) and neutral kaon mixing param- 
eter e for different representative values of the Higgs mass {Mh = oo, red triangle; Mh = 75 
TeV, blue square; Mh = 20 TeV, large green dots). Reprinted from Nucl. Phys. B., 802, Y. 
Zhang, H. An, X. Ji, R. N. Mohapatra, "General CP violation in minimal leftright symmetric 
model and constraints on the right-handed scale", p. 247., Copyright (2008), with permission 
from Elsevier . 

6.3 Additional Spacetime Dimensions 

The study of EDMs in BSM scenarios involving extra spacetime dimensions, such as the Randall- 
Sundrum (RS) paradigm for warped extra dimensions or flat but orbifolded extra dimensions, is con- 
siderably less advanced than in the case of SUSY or extended gauge symmetries. To our knowledge, 
EDM estimates have been largely confined to the use of NDA, coupled with an analysis of the flavor 
and CP structure associated with a given implementation of the RS paradigm. 

For concreteness, we focus on the scenario RSI, wherein SM fields may propagate in the "bulk" 
of the fifth dimension between two branes: the TeV brane and the Planck brane. The Higgs field is 
localized at the former while gravity lives at the latter, ensuring that the natural scales for the EW and 
gravitational interactions are the weak and Planck scales, respectively. The dependence on all mass 
scales in the fifth dimension follows from the "warping" associated with the dependence of the metric 
on the fifth dimensional co-ordinate z as 



1 



(6.230) 



(kz) 
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where denotes the usual four dimensional co-ordinate vector and k is a warping factor. 

Flavor structure arises from the 2;-dependence of the fermion wavefunctions rather than from the 
values of the Yukawa couplings for the 5-dimensional theory. Assuming the latter to be "anarchical", 
the observed fermion mass hierarchy arises when the light fermions are localized near the Planck brane 
and the top quark near the TeV brane. Since the Higgs is localized near the TeV brane, its vev gives a 
significantly larger mass to the top quarks than to the light fermions. The Kaluza-Klein (KK) modes 
for the light fermions are also localized near the TeV brane, generating a leading order suppression of 
flavor changing neutral currents (FCNCs). 

An early concrete application of this scenario to the flavor and CP problems was carried out in 
Ref. |159] . The corresponding five-dimensional Lagrangian contains two components of interest: 

(a) the bulk Lagrangian 

Aermion = {^F^^Z^mV^ + kCQUQud)iQud)} (6.231) 

where all of the SM fields propagate in the warped extra dimension (denoted by a co-ordinate z), 
where Dm is the five dimensional covariant derivative, and where Cg^d are 3x3 Hermitian matrices 
that determine the 5-D masses. 

(b) the 5-D Yukawa interaction: 

-Cbrane = hS{z - Z^)\l'^^Q{u, d) (6.232) 

where A^^ are the Yukawa matrices, h is the Higgs field, and zq indicates the location of the TeV-brane. 

Carrying out the Kaluza-Klein (KK) reduction of the 5D theory to an effective 4D theory on the 
TeV brane yields the SM fields (zero modes) and their KK partners. The quark zero mode masses 
are then given by niq ~ vFqX^^^Fu^d, where FQ^u,d are the values of the quark wavefunctions on the 
TeV brane . For purposes of the present discussion, the specific values of the FQ^u,d are not essential. 
However, due to the different profiles for the light fermion zero mode and KK modes, couplings between 
the two go as A^^Fg. Since these couplings are not aligned with the quark masses, non-trivial flavor 
and CPV contributions may be generated at one-loop order. Representative diagrams that generate the 
quark EDMs are shown in Fig. |9| In each case, an odd number of zero mode - KK Yukawa interactions 
is needed to obtain the chiral structure associated with the EDM. The gauge loops contain only one 
insertion while the Higgs loops contain three. In the former case for a down quark-gluon loop, one has 

drffeluon, KK] ~ kv (^DlFgXfFdDn) ~ [diag(mrf, m„ rrib)]^^ = , (6.233) 

where Dlr rotate the left- and right-handed down quarks between the flavor and mass bases. In 
contrast, the Higgs loop results are not aligned with the light quark mass matrix and, as shown in 
Ref. ^159j . lead to a non- vanishing EDM contribution: 

(6.234) 

Taking the phases that enter this expression to be maximal, the authors of Ref. |159j arrive at the NDA 
estimate 

.„1H.^, KK, . I (^) i^-^) . ,.-"e M X (^)^ (IM)^ ,,.35, 

An earlier analysis by the authors of Ref. |160] that included the contribution from the CEDM found 
a considerably smaller sensitivity to the CPV parameters. In that study, the additional suppression 
results from a tiny coupling between the first and third generation as well as constraints from the CPV 
parameter that enters the neutral kaon system. 



rfdpiggs, KK] ~ 2k^v 
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In addition to the loop contribution, one expects contributions from operators that hve on the 
TeV brane. This term is UV-sensitive and, thus, depends on the cutoff of the effective theory on the 
TeV brane, obtained from the cutoff at the Planck scale by warping: A ~ As^ exp(— vrfcrc), where 
Trkvc ~ Mpianck/TeV to solve the hierarchy problem. From NDA one anticipates 

d„[brane] ~ eC^ (^) ~ (IQ-^^ e fm) x {^^^ ^^-^^e) 

The value of A depends on the specific realization, depending on whether the Higgs is placed on the 
TeV brane, in the bulk, or in the bulk but localized near the TeV brane. For sin0cpv ~ ^^(1); the 
resulting contribution to dn can be comparable to the present experimental limit or comparable to the 



considerably larger loop contribution in Eq. (6.235). 



One should bear in mind that the foregoing results are obtained using NDA and that loop com- 
putations in extra-dimensional scenarios are subject to theoretical ambiguities. Nonetheless, one thus 
finds a situation similar to that in SUSY: current EDM limits imply that either the CPV phases are 
suppressed or that the KK mass scale lies well above the TeV scale. A variety of solutions to the RS 
CP problem have been proposed. In Ref. |161j a variant of RSI was analyzed under the assumption of 
5D minimal flavor violation, leading to the vanishing of EDMs at one-loop order. Ref. [33] considered 
an RSI scenario with spontaneous CPV, where the source of CPV was geometrically sequestered from 
the TeV brane by placing it in the bulk. The model provides a natural suppression of 6, while the 
dimension-six EDM operators first appear at two-loop order. 

These studies notwithstanding, it is evident that there exists considerable room for further work on 
EDMs in extra dimensional models. To our knowledge, no computations of the CEDM, three-gluon, 
or dimension-six four-fermion operators has appeared in the literature. In particular, the limits on 
dAC^^Hg) may imply more severe constraints on RS CPV than have been obtained in these earlier 
studies, given the long-range vr-exchange contributions to the nuclear Schiff moments sourced by the 
CEDM operators. Even with the ambiguities associated with loop computations in c? > 4 dimensions 
and with cutoff-dependent TeV brane operators, a study of these additional CPV effects would be both 
interesting and potentially significant. 



7 Outlook 

In the context of fundamental symmetry tests during the LHC era, searches for the permanent electric 
dipole moments of atoms, molecules, nucleons and nuclei will provide one of the most powerful probes 
of both BSM physics as well as the remaining, as-yet unseen source of SM CPV - the QCD ^-term. 
Improvements in experimental sensitivity are poised to improve by as much as two orders of magnitude 
in the near term and possibly further on a longer time scale. The observation of a non-zero EDM 
would constitute a major discovery, pointing to a non- vanishing 6 parameter and/or a new source 
of CPV associated with new fundamental interactions. Conversely, the non-observation of EDMs at 
the anticipated sensitivity levels would tighten the already stringent those on 6 as well as those on a 
variety of BSM scenarios. Either way, the implications for the fundamental laws of nature and their 
consequences for the cosmic baryon asymmetry cannot be overstated. 

Theoretically, the challenge is to provide the most robust framework for interpreting the results of 
EDM searches and delineating their implications. Doing so entails analyzing physics associated with a 
variety of energy scales, ranging from the short-distance physics of CPV at the elementary particle level 
to the longer-distance physics at the hadronic, nuclear, and atomic/molecular scales. In this review, we 
have relied on effective CPV operators at mass dimensions four and six as a bridge between the physics 
of these various scales. While the use of effective operators is not applicable in all circumstances, such 
as those in which BSM CPV involves new light degrees of freedom, it nevertheless provides a broadly 
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Figure 9: Representative contributions to the d-quark EDM and CEDM in the Randall- 
Sundrum scenario. 



apphcable and model-independent context for the interpretation of EDM experiments. Limiting our 
consideration photons, gluons, and first generation fermions, we encounter thirteen presently undeter- 
mined CPV parameters at d = 4 and d = 6: 6, the fermion EDMs, quark CEDMs, three-gluon operator, 
and several four-fermion operators. The task for theory, then, is to delineate how these operators may 
be generated by BSM physics above the weak scale, how they evolve to the hadronic scale, and how 
they generate the appropriate hadronic, nuclear, atomic and molecular matrix elements that ultimately 
give rise to EDMs in these systems. 

From our review of this theoretical effort, several features emerge: 

(i) The EDMs of paramagnetic atoms and molecules are dominated by two quantities: the electron 
EDM and one combination of semileptonic, CPV four-fermion operators characterized by the 
Wilson coefficient Im Ciq \ Moreover, when characterizing the former in terms of the appropriate 
dimensionless parameter 6e, one finds that the EDMs of these systems are an order of magnitude 
more sensitive to Im Ceq^ than to Se- The level of theoretical atomic/molecular theory uncertainty 
in either sensitivity is roughly 10% or better. For Im Ceq \ the associated hadronic matrix element 
of the isoscalar scalar density is under reasonable control, given that it can be obtained from the 
pion-nucleon a-term and the average light quark mass. 

(ii) Diamagnetic atom EDMs are most sensitive to the nuclear Schiff moment, individual nucleon 
EDMs, and the semileptonic four-fermion tensor operator with Wilson coefficient Im C\'. 
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(iii) Neither system provides a particularly sensitive handle on Im Ceq , given the relative suppression 
of the associated isovector scalar or isoscalar pseudoscalar nuclear matrix elements. 

(iv) There exists considerable room for refinement in computing the nucleon matrix of underlying CPV 
operators as well as nuclear Schiff moments. While the sensitivity of g^^ to 6 and the dependence 
of (i AT on the quark EDMs is now known fairly reliably, the uncertainties associated with matrix 
elements of the quark CEDMs, three-gluon operator, and four-quark operators are large. Similar 
statements apply to the dependence of c? at on ^ as well as to the sensitivity of the nuclear Schiff 
moments to the g^^ . Perhaps, one of the primary challenges facing is now to achieve a more 
reliable set of hadronic sensitivities. 



(v) Searches for the EDMs of diamagnetic atoms and nucleons alone is unlikely to disentangle the 
effects of the quark CEDMs, three-gluon operator, and four-quark operators. However, the pos- 
sibility of searching for EDMs of the proton, deuteron, triton and helion in storage rings would 
offer additional handles on these underlying sources of hadronic CPV based on their chiral trans- 
formation properties. 

(vi) Any global analysis of EDM search results, whether performed at the model-independent level 
of the effective operators within a given BSM scenario such as supersymmetry, should take into 
account the rather sizeable theoretical hadronic and nuclear uncertainties associated with the 
sources of hadronic CPV. 



Setting aside the aforementioned uncertainties, current EDM null results imply that any new CPV 
lies at the TeV scale or that CPV phases are O{10~^) in magnitude or smaller. The next generation of 
searches will push these sensitivities to A > 10 TeV or equivalently | sin0cpv| ~ C(10~^), putting their 
reach well beyond that of the LHC. Should the LHC have observed only the SM Higgs boson (a major 
discovery in its own right) by the end of this decade, then EDM searches will provide one of the most 
effective tools for probing the next piece of terrain in the high energy desert. Given these prospects, 
addressing the open theoretical challenges summarized above becomes all the more important. We hope 
that our discussion of the rich array of physics associated with EDMs will spur new efforts to take on 
these challenges. 
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A Scalar and pseudoscalar form factors: heavy quark contri- 
butions 

The g'^ for the heavy flavors can be obtained by using the trace anomaly for the energy momentum 
tensor and nucleon mass and by integrating out the heavy quarks \162\ 1163] . We first have 



^m^Pn = {N\ \N) = ^ {N\ m,qq \N) + ^ (iV| Tr {G.^G^ \N) , (A.237) 



q=u,d,. 



2a. 



where 



j3(as) = -9a^/27r (A.238) 
and where we have used the result that integrating out each heavy quark leads to the replacement 

2 a. 



{N\mQQQ\N) = -^- ^ {N\TTiG,.Gn\N) • 



(A.239) 



Using Eqs. ( |3.92[3.100| ) we may solve for (A^| QQ \N) and, thus, g^- 

9a 



47r 



{N\ Tr {G^^n \N) = [m^ - (m^v)^ - 220«:,] ^n^n , 



(A.240) 



so that from Eq. ( |A.239D we obtain Eq. ( [3l00| . 

For the gp of heavy flavors, we follow |162l [T2| [62] and exploit the U(l)yi anomaly. Letting 

<//.5 = ^^A'^s? + Yl ^^m75Q (A.241) 



q=u,o 



Q=s,c,b,t 



and 



we have 



{N\J^,\N)=gf {N\N^,^,N\N) 



{N\ d^J^, \N) = 2m^(?f (ATI m^,N \N) = 2 (^1 "^.^^75? \N) 



(A.242) 



(A.243) 



q=u,d 



+2 Y {N\mQQt^,Q\N)+Q{N\^Tr(G,,G^'')\N) 

Q=s,c,b,t 



Using 



{N\ rriQQi^.Q \N) = (A^| Tr (G,,G^^^ \N) 
and the expressions for matrix elements of ui'-f^u ± di'y^d in terms of g^g'^^ we obtain 



(A.244) 



{N\mQQt^,Q\N) = 

{N\ 



{0)mN , (o)2m 

9 a +9p 

ruQ 



1 (1) Anig . . 
4 niQ 



Q{0) Q(l) 

gp + gp 7-3 



(A.245) 
(A.246) 



with the result in Eq. (3.104). 
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B CPV parameter sensitivities: a compilation 
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Cn 


8.2 X 10-9 


(4 - 12) X 10-9 










^d 


-33 X 10-9 


-(16-50) X 10-9 
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Pi 


2 X 10-^ 


(0.2-40) X 10-^ 


7(^) 


2 X 10-6 


(1 - 10) X 10-6 


Im C<^ud 




3 X 10-^ 


(1 - 10) X 10-^ 


(pud 

7(1) 


1 X 10-6 


(5 - 150) X 10-^ 


Tm r^^'^) 

^quqd 


pquqd 


40 X 10-^ 


(10 - 80) X 10-^ 


quqd 

7(i) 


2 X 10-6 


(1 - 10) X 10-6 




ft?' 


12.7 


11-14.5 












0.9 


0.6-1.2 









Table 6: Best values and reasonable ranges for hadronic matrix elements of CPV operators. First 
column indicates the coefficient of the operator in the CPV Lagrangian, while second column indicates 

the hadronic matrix element (sensitivity coefficient) governing its manifestation to the neutron EDM. 
Third and fourth columns give the best values and reasonable ranges for these hadronic coefficients. 
Firth to seventh columns give corresponding result for contributions to TVPV ttNN couplings. " Units 
are e fm for all but the and p^. ^ We do not list entries for (7(^)^, (^(i), tKj.)) ^ they are suppressed by 
ck/tt with respect to (7^7, ti^(i), fl{i)) ■ The W(o,i) are in units of fm-^. 
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LEG 


do 


di 




gir efm 


(F^COefm 


(F3C2)efm 


dn 


1 


-1 


- 


- 


- 


- 


dp 


1 


1 


- 


- 


- 


- 


dd 


2 





-0.0002 + 0.07/3i 


0.2 


- 


- 


dh 


0.83 


-0.93 


0.1 


0.2 


-0.01 


0.02 


dt 


0.85 


0.95 


-0.1 


0.2 


0.01 


-0.02 



Table 7: Dependence of the EDMs of the neutron, proton, deuteron, hehon, and triton on the six 
relevant TVPV low-energy constants at leading order. A "-" denotes that the low-energy constant 
does not contribute in a model-independent way to the EDM at this order. For the potential-model 
dependence and other uncertainties in the results, see text. The P-, T-even isospin-breaking pion- 
nucleon coupling /3i is not well known, /3i = (0 ± 9) ■ 10"^ [ ^1113] . (Adapted from Ref. |103j .) 



Source 


^^-term CEDM quark EDM chiral-invariant left-right 


Ay^dn/e 
dp/ dn 
dd/ dn 
dh/ dn 
dt/dh 


0(§ImC,) 0(^ImC(i'3) 
Oil) Oil) Oil) Oil) Oil) 
Oil) 0{AI/Q^) Oil) Oil) 0{Al/Q^) 
0{Al/Q') 0{Al/Q') Oil) Oil) 0{Al/Q') 
Oil) Oil) Oil) Oil) Oil) 



Table 8: Expected orders of magnitude for the neutron EDM (in units of e/A^), and for the EDM 
ratios proton to neutron, deuteron to neutron, helion to neutron and triton to helion, for the ^-term 
and dimension-six sources. For chiral-invariant sources, ImC^ stands for Cq and ImC^^'^j. Q represents 
the low-energy scales F^^, m^r, and A/m/v~B, with B the binding energy. (Adapted from Ref. |103j .) 



Ref. 


Method 


ao 




02 


m 


Schematic 


0.087 


0.087 


0.174 


m 


Phenomenological RPA 


0.00004 


0.055 


0.009 


m 


Skyrme QRPA 


0.002 - 0.010 


0.057 - 0.090 


0.011 - 0.025 


[125j 


Odd-A Skyrme mean-field theory 


0.009 - 0.041 


-0.027 - +0.005 


0.009 - 0.024 



Table 9: The coefficients aj in ^^^Hg from a variety of nuclear-structure calculations. 



Ref. 


Method 


Go CLl 0'2 


[129j 

m 


Octupole-deformed Wood-Saxon 
Odd-A Skyrme mean-field theory 


-18.6 18.6 -37.2 
-1.0 --4.7 6.0-21.5 -3.9- -11.0 



Table 10: The coefficients in ^^^Ra from two nuclear-structure calculations. 
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Ref. 


Nucl. 


Method 


«o 


ai 


a2 


m 




Schematic 


-0.11 


-0.11 


-0.22 


m 




Phenomenological RPA 


-0.008 


-0.006 


-0.009 


m 




Phenomenological RPA 


0.019 


-0.061 


0.053 


[125J 




Odd-A Skyrme mean-field th. 


0.034 - 0.042 


-0.0004 - -0.018 


0.064 - .071 


m 




Phenomenological RPA 


0.012 


0.021 


0.016 


[129] 


223-D „ 

Ra 


Octupole-def. Wood-Saxon 


-lb 


25 


- 50 


[129J 




Octupole-def. Wood-Saxon 


-62 


62 


- 100 


wm 


223p^ 


Octupole-def. Wood-Saxon 


-31 


31 


- 62 


\m 


219Fr 


octupole-quadrupole vibr. 


-0.02 


-0.02 


-0.04 



Table 11: The coefficients in some other nuclei of interest, from several nuclear-structure calculations. 



Nucl. 


Best value 
ao cti (22 


Range 

CtQ (^1 0,2 


199gg 

i29Xe 

225Ra 


0.01 ± 0.02 0.02 
-0.008 -0.006 -0.009 
-1.5 6.0 -4.0 


0.005 - 0.05 -0.03 - +0.09 0.01 - 0.06 
-0.005 - -0.05 -0.003 - -0.05 -0.005 - -0.1 

-1 - -6 4 — 24 -3 - -15 



Table 12: Best values and ranges for the coefficients in three nuclei used in experiments. 



Atom 


P\ 


P7 


<1 


Pp 


Pn 


/«5 






e fm 


10-*^ e fm 


10-4 


10-4 


10-4 fm-2 


205 rjj 


-573(20) 


0.65(0.02) 


189(7) 








133CS 


123(4) 


-0.14(0.005) 


-41(1.3) 








85Rb 


25.7(0.8) 


-0.03(0.0009) 


-8.5(0.3) 








210p^ 


903(45) 


-1.02(0.05) 


-298(15) 








199 Hg 


0.01 


-1.13 X 10-5 


-3.3 X 10-3 


-0.56 


-5.3 


2.8 


Molecule 


Hz/(e fm) 


Hz 


kHz 








YbF 


-(1.1 + 0.1) X 10^2 


(1.2 + 0.1) X 10^ 


(3.6 + 0.3) 








ThO 


-(4.6 + 0.4) X 10^2 


(5.2 + 0.5) X 10^ 


(15.2 + 1.4) 









Table 13: Dependence of atomic and molecular EDMs on EDMs of the electron, proton and neutron 
and on the Schiff moment. 
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Table 14: Summary of the combined bounds at 95% c.l. on three phases (02, 03, 4>u,d) for tan/3 = 3, 
60 and first generation sfermion masses {Mll^rrJi = 200, 500, and 1000 GeV, obtained using current 
experimental limits on the neutron. Thallium, and Mercury EDMs [TM] .With kind permission from 
Springer Science+Business Media: J. High Energy Physics, "A comprehensive analysis of electric dipole 
moment constraints on CP-violating phases in the MSSM", 08, 2010, p. 062, Y. Li, S. Profumo, M.J. 
Ramsey-Musolf , Table 4. 



tan/3 


3 


60 


{Mll,rr)i 


200 GeV 


500 GeV 


1000 GeV 


200 GeV 


500 GeV 


1000 GeV 


I02| 




< 2.1 X 10-3 


< 5.0 X 10-3 


< 1.5 X 10-2 


< 9.3 X 10-5 


< 2.5 X 10-^ 


< 6.9 X 10-^ 


I03| 




< 2.8 X 10-3 


< 9.7 X 10-3 


< 2.8 X 10-2 


< 3.1 X 10-^ 


< 4.2 X 10-^ 


< 1.5 X 10-3 


\<Pu,d 




< 1.8 X 10-2 


< 6.0 X 10-2 


< 0.17 


< 1.7 X 10-2 


< 5.6 X 10-2 


< 0.21 



CPV Parameter 


Coefficient 


Method 


Value 


Remarks 


e 


an 


ChPT 


~ 0.002 e fm 


See Eq. ( 


3.75 


) 


e 




Lattice QCD[79] 


-0.040(28) e-fm 


= 0.53 GeV 


e 


(y.p 


Lattice QCD[79] 


0.072(49) e-fm 


= 0.53 GeV 


e 




Lattice QCD[SU] 


-0.049(5) e-fm 


^ 0.61 GeV 


e 


Op 


Lattice QCD[8n] 


0.080(10) e-fm 


^ 0.61 GeV 


e 


an 


QCD Sum Rules [82 W\ 


(0.0025 ± 0.0013) e-fm 


A from QCD SR 


e 




QCD Sum Rules[H2] 


(0.0004inS^^) e-fm 


A from lattice 


e 


A(o) 


ChPT 


~m2/A^F^ -0.08 
0.017 ±0.005 


See Eq. ( 
See Eq. ( 


3.71 
3.70 


) 
) 


e 


A(i) 


ChPT 




See Eq. ( 


3.73 


) 



Table 15: Dependence of hadronic quantities on 6. 
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CPV Parameter 


Coefficient 


Method 


Value 


Remarks 


Ini Cj„r< 


QuG 


ChPT 


~ e/f ~ 8 X 10~^ e fm 




Trn Cj,,r< 


QuG 


OCD SRI85 


-fO 6 ± 3) X 10-3 e fm 


PO assumed 


Im Cm 

^^^^ ^ a{j 


QdG 


QCD SR[85J 


-(1.2 ± 0.6) X 10-3 e-fm 


PQ assumed 


Trn Cj,,r< 


ouG 
h'n 


OCD SRE2] 


-(0 20+nn^) X 10-3 e-fm 


PQ assumed 


Im C^n 


QdG 


QCD SR[82] 


-(0.40+0-33,) X 10-3 e-fm 


PQ assumed 


ImC r- 

±111 \-yq(J 


QuG 


OM/NDA 


~ 1 X 1 0~^ p fm 

X /N X V7 l^' ±±±± 


iurluHps hC /-< 

llH^l LHJ-V^O ^^CtKj 


Im C„r^ 


QdG 


QM/NDA 


4 X 10-"^ e fm 


includes KqQ 


d„ 


Oat 


ChPT 


~ -0.7 






oil 


OCD SRf85l 


55 ± 2 8 


PO assumed 


dn 


pii 


QCD SR[85] 


1.1 ± 0.55 


PQ assumed 


d 




OM/NDA 


~ —0.09 






HN 


QM/NDA 


~ 0.36 


iiicludps K^n 


L 
"<? 




ChPT 


~ 5 X 10-^ e fm 








OCD SRI85I 


-(0 9 ± 5) X 10-*^ e fm 


PO assumed 




eCi 


QCD SR[85] 


(-3.6 ± 1.8) X 10-*^ e fm 


PQ assumed 






QM/NDA 


~ 0.2 X 10-s e fm 


includes K„c 


Sd 




QM/NDA 


0.8 X 10-*^ e fm 


includes Kqc 



Table 16: Dependence of nucleon EDM on quark CEDMs expressed in terms of the quantities ImCgc, 
dq, or 5q. 



CPV Parameter 


Coefficient 


Method 


Value 


Remarks 


lm{CuG±CdG) 




Chiral/NDA 


~ 0.03 




Im {CuG + Cdc) 




QCD SRIM] 


(-1.7 o 0.6) X 10-2 


PQ assumed 


Im {CuG - CdG) 




QCD SR[H6] 


-(2.3^^1) X 10-2 


PQ assumed 


dq 


^{0,1) 


Chiral/NDA 


~ -26 




du + dd 


W(0) 


QCD SR[86j 


(-5 ^ 15) fm-^ 


PQ assumed 


du dd 




QCD SR[86] 


20;to fm-i 


PQ assumed 




mi) 


Chiral/NDA 


~ 1.9 X 10-6 




^u + ^d 


^(0) 


QCD SR^ 


(-3.5 ^ 1.2) X 10-^ 


PQ assumed 


5u - 5d 


^(1) 


QCD SRp] 


(-4.6^^:^) X 10-^ 


PQ assumed 



Table 17: Dependence of TVPV ttNN coupling on quark CEDMs expressed in terms of the quantities 
CqG, dq, or S. 
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CPV Parameter 


Coefficient 


Method 


Value 


Remarks 


C 


In 


Chiral/NDA 


~ e/v ^ 8 X IQ-"^ e fm 




c 


ou-y 


QCD Sum Rules[Il] 


(0.4 ±0.2) X 10-3 e fm 




Cd-y 




QCD Sum Rules[Tl] 


-(1.6 ±0.8) X 10-3 e fm 




Cu-y 


ouy 


QCD Sum Rules[82] 


(0.13+°:J°) X 10-3 e fm 




Cdj 


Pf 


QCD Sum Rules[H2] 


-(0.531°;^^) X 10-3 e fm 




Cwy 


ouy 


Quark Model 


0.4 X 10-3 e fm 




Cd-y 


Pic 


Quark Model 


— 1.5 X 10-3 e fm 




c 


ou-y 

yn 


PQM[T2l|8H] 


0.6 X 10-3 e fm 




Cd-y 


Pi? 


PQM[I21EH] 


-0.8 X 10-3 e fm 




dg 


Pn 


Chiral/NDA 


~ -0.7 




du 


Pi 


QCD Sum Rules[Tl] 


(-0.35 ±0.17) 




dd 


pi 


QCD Sum Rules[Il] 


(1.4 ±0.7) 




du 


Pi 


QCD Sum Rules[H2] 


(-o.iilSD 




dd 


pi 


QCD Sum Rules[S2] 


(0.47l°:f) 




du 


Pi 


Quark Model 


-1/3 




dd 


pi 


Quark Model 


4/3 




du 


pi 


PQM[I21IHH] 


(Am)„ = -0.508 




dd 


pi 


PQM[T2llH8] 


(Arf)„ = 0.746 




ds 


pi 


PQMinilHH] 


(A.)„ = -0.226 






<l 


Chiral/NDA 


~ 5 X 10-8 e fm 




K 




QCD Sum RulesHg 


(0.8 ± 0.3) X 10-8 e fm 




5d 


<i 


QCD Sum Rules[Il] 


(-3.2 ±2.3) X 10-8 e fm 




5u 




QCD Sum Rules[H2] 


(0.271°:?) X 10-8 e fm 




Sd 


< 


QCD Sum Rules[82] 


(-l.ll^l) X 10-8 efm 




Su 




Quark Model 


0.8 X 10-8 e fm 






eCn 


Quark Model 


-3.2 X 10-8 e fm 




5u 


eC 


PQMinilHH] 


1.2 X 10-8 e fm 




5s 


eC 


PQM[T2l|88] 


-1.7 X 10-8 e fm 





Table 18: Dependence nucleon EDM on quark EDM expressed in terms of the Wilson coefficients 
Cq^(A^), individual quark EDMs dg{A^), or dimensionless quantities 6g{A^). 
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UrV rarameter 


Uoemcient 


Method 


Value 


illi ^ Q 




PViirnl /ND A 

v^iiii ai / IN 


~ 4(1 V 1 p fm 






QCD SR[87] 


2.0 X 10-^ e fm 




7(i) 


Chiral/NDA 


~ 2 X 10-6 



Table 19: Dependence of hadronic quantities on Weinberg three-gluon operator Wilson coefficient. 



CPV Parameter 


Coefficient 


Method 


Value 


Remarks 


^quqd 

Im C^ud 


^^ud 
P^ud 


Chiral/NDA 
Chiral/NDA 
Saturation /ChPT 


~ 40 X 10-^ e fm 
~ 3 X 10-s e fm 
1.3 X 10"^° e fm 


NLO 


Im C^ud 
\Ya.C^ud 


^(0,1) 

^(1) 

Lpud 
^(1) 


Chiral/NDA 
Chiral/NDA 
Saturation 


~ 2 X lO"*' 

~ 10-6 
3.3 X 10-5 





Table 20: Dependence of hadronic quantities on CPV four-quark operators. The Saturation result for 
Z?^"^ has been obtained by first computing g^"* and then employing the NLO result from ChPT given in 
Eqs. (|3l07pl08|. 



Form Factor 


Method 


Value 


Remarks 


sf 

9f 
9? 


Lattice QCD 
Lattice QCD 
Lattice QCD 

Chiral 

Chiral 
Lattice QCD 
Lattice QCD 


6.3 ±0.8 
0.45 ±0.15 
0.4 (2) 

0.53(18) 


Eq. (3.92) and [53l [M] 
Eq. (3.92) and [58] 
Isospin and [63] 

Isospin and [M] 



Table 21: Form factors entering semileptonic CPV interactions. 
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System 


Present 90 % C.L. 
Limit (e fm)" 


Sensitivity Goal^ 


Group 


SM CKM (e fm)^ 


Cs 


1.2 X 10-1° 




|164j 


~ 10-23 


Tl 


9.5 X 10-12 






~ 10-22 




10.5 X 10-15 




m 


~ 10-19 


ThO'^ 


- 


10-15 10-1^ 






n 


2.7 X 10-13 




[166] 


1.6 X 10-1^ ^ 1.4 X 10-20 


n 




(1 -3) X 10-1^ 


CryoEDM 




n 




4 X 10-15 


nEDM/SNS 




n 




5 X 10-1^ 


nEDM/PSI 




n 




5 X 10-15 


n2EDM/PSI 




n 




2 X 10-15 


nedm/FRM-II Munich 




n 




10-14 _ 10-15 


TRIUMF 




P 




10-16 


srEDM 




199gg 


2.6 X 10-16 


(2.6-5) X 10-1^ 


m 




225Ra 




(10 - 100) X 10-15 


Argonne 








1.3 X 10-1^ 


TRIUMF 








2 X 10-15 


FRIB 




223Xe 


5.5 X 10-1^ 




[168] 





Table 22: Present EDM limits and sensitivity goals for the paramagnetic atoms and molecules (first 
group); nucleons {n,p) (second group); and diamagnetic atoms (third group). A limit on the electron 
EDM of 10.5 X 10"i5 e fm (90% C.L.) has been derived from the most recent YbF experiment |147j 
assuming it is the only source of the molecular EDM (see Section [s]). Also listed are the expected 
magnitudes of the SM "background" due to the phase in the CKM matrix. " We thank T. Chupp 
for providing the 90% C.L. limits from existing searches. All sensitivity goals are self-reported by 
members of the given collaboration. We do not quote SM CKM predictions for diamagnetic atoms, 
due to the incorrect chiral implementation of chiral symmetry in Ref . |114j as pointed out in Ref . |169] . 
Molecular sensitivity expressed in terms of limit on de rather than on (Ia- 
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Atom 


l(0) 


^.(0) 

fVp 


1.(0) 

fVij-i 




10-^ e fm 


10-^°e fm 


10-^ e fm 


205 


-7.0 ±0.3 


-1.5 


-0.5 




-0.78 ±0.2 


-2.2 


-0.92 




-0.110 ±0.003 






210pj. 


-10.9 ± 1.7 






199Hg 


-8.1 X 10-^ 


6 


4 


Molecule 


kHz 


kHz 


kHz 


YbF 


-92 ±9 






ThO 


-561 ± 56 







Table 23: Dependence of atomic and molecular EDMs on T- and P-odd semileptonic interactions. 
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